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ABSTRACT 


During  the  three  year  period  under  the  grant,  01  Feb  1985 
to  31  Jan  1988,  telescopic  observations  were  carried  out  on 
selected  infrared  sources,  theoretical  investigations  were  pur¬ 
sued,  five  papers  were  published,  a  sixth  had  been  submitted, 
and  work  was  going  forward  at  the  end  of  the  grant  period  on 
the  investigation  of  the  infrared  nature  of  planetary  nebula 


nuclei 


I.  SYNOPSIS 

Start  date  on  the  grant  was  01  Feb  1985,  but  funding  was  not 
made  available  by  AFOSR  until  after  15  April  1985.  Computer  equip¬ 
ment  purchased  with  grant  funds  was  acquired  and  placed  in  service 
during  the  summer  of  1985.  At  the  same  time  a  student  assistant, 
David  Mikelson,  was  hired  to  assist  with  telescopic  observations, 
computer  programming,  and  data  reduction. 

The  announcement  in  the  fall  of  1985  by  AFGL  that  the  LAIRTS 
program  would  not  be  funded  had  a  major  impact  on  the  scientific 
program  to  be  carried  out  under  this  grant.  LAIRTS  science  plan¬ 
ning  was  discontinued,  and  major  emphhsis  was  shifted  to  ground 
based  infrared  observing  and  analytical  studies. 

In  early  1986  David  Mikelson  resigned  from  his  position  as 
assistant,  and  for  the  remainder  of  the  contract  period  I  was 
unable  to  recruit  a  qualified  replacement  for  him  from  among  the 
EMC  student  body.  Eastern  Montana  College  is  a  poorly  funded 
state  institution,  which  in  some  respects  has  not  evolved  into 
a  true  four-year  college  with  a  balanced  Arts  and  Science  curri¬ 
culum,  with  the  result  that  we  have  very  few  science  oriented 
students  on  campus. 

Purchase  of  IRAS  data  was  initially  postponed  when  it  was 
learned  that  the  data  set  would  cost  approximately  $4000,  as 
compared  to  $200  estimated  in  the  original  cost  proposal.  In 
early  1986  unexpended  funds  from  the  first  year  of  the  grant 
became  available,  and  permission  to  use  a  portion  of  these  funds 
for  purchase  of  the  IRAS  data  was  requested  from  AFOSR.  After 
a  full  six-month  delay,  this  request  was  granted.  By  then  it 
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had  become  apparent  that  qualified  student  help  was  not  available, 
so  I  put  off  ordering  the  IRAS  data  until  I  had  my  student  help 
problem  solved.  This  never  happened,  so  the  IRAS  database  was 
never  purchased. 

Despite  the  delays  and  problems  mentioned  above,  research 
under  the  grant  went  forward  in  terms  of  continuing  scientific 
observations  and  the  publication  or  preparation  for  publication 
of  research  results. 

As  a  result  of  the  lack  of  qualified  student  assistants  on 

this  campus,  a  paucity  of  meaningful  College  support  for  scien- 

« 

tific  endeavor,  and  a  recent  administrative  decision  to  de-empha- 
size  the  physical  sciences  even  further,  I  have  regretfully  decid¬ 
ed  against  proposing  further  AFOSR  grant  support  for  my  research. 


RESEARCH  STATUS 

A.  Objectives 

The  scientific  objectives  of  this  research  program,  as 
stated  in  the  grant  proposal,  were: 

(1)  to  determine  the  nature  of  infrared  sources  observed. 

(2)  to  establish  the  statistical  distribution  of  infrared 

-20 

sources  in  space  to  a  sensitivity  limit  of  10 

_2 

W  cm  and  a  spatial  resolution  of  4  arcseconds. 

B.  Revised  Objectives 

Due  to  cancellation  of  the  LAIRTS  program  by  the  Air 
Force,  the  follo&iiig;  revised  objectives,  for  this  re¬ 
search  progranKwere  formulated: 


formulated: 


-XCO  to  determine  the  nature  of  infrared  sources  observed. 
/  / 


t 

* 

* 

l 


f 
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(2) '"‘to  carry  out  theoretical  investigations  on  the  nature 

and  evolution  of  selected  infrared  sources.  4 

/ 

^£7)  to  establish  the  statistical  distribution  of  infra¬ 
sources  in  space  to  a  sensitivity  limit  and  spatial 
resolution  consistent  with  available  observing 
equipment  and  techniques.  ^ 

C .  Research  Progress  / 

During  the  grant  period,  01  Feb  1985  to  31  Jan  1988, 
progress  towards^the  research  objectives  was  made,  as 
■jfollow^; 

(1)  ^Research  on  infrared  sources  in  the  W40  complex 
was  published, in  the  Astro physical  Journal,  a 
copy  of  the  paper  is  included  in  the  appendix. 

~~~ — (-2-}-}  Results  of  radio  and  infrared  observations  of 

OH/IR  stars  in  the  Galaxy  were  published,  in  the 
Astronhysical  Journal.  Acopy  of  this  paper  is 
included  in  the  appendix. 

.j  _ ( 3 ) j  A  theoretical  study  of  the  possibilities  for  obser¬ 

vation  of  galaxies  in  the  process  of  formation  in 
the  early  universe  was  published  in  the  Proceedings 

of  the  Montana  Academy  of  Science.  A  copy  of  this 

. 

paper  is  included  in  the  appendix. 

(4)  In  collaboration  with  a  number  of  colleagues  from 
Wyoming^  the  results  of  a  series  of  observations 
made  on  Nova  Vulpeculae  1984  fr'2  were  published  in 
the  Astronhysical  Journal  (Letters).  A  copy  of 


this  paper  is  included  in  the  appendix 


w 


(5)  Preliminary  results  on  our  studies  of  the  infrared 
properties  of  the  nuclei  of  planetary  nebulae  were 
presented,  in  Mexico  City  at  the  IAU  Symposium  #131 
on  planetary  nebulae,  in  October  1987.  An  abstract 
of  this  paper  is  being  published  in  the  Symposium 
Proceedings,  and  a  copy  is  included  in  the  appendix. 

(6) ^A  theoretical  model  has  been  constructed  for  the 

transfer  of  cosmic  ray  energy  to  emission  nebulae  , 

and  has  been  sent  to  the  Astro physical  Journal 

for  publication,  and  the  paper  has  been  returned, 

« 

with  reviewer's  comments,  for  revision.  A  copy  of 
the  original  draft  of  this  paper  is  included  in  the 
appendix. 
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INFRARED  SOURCES  AND  EXCITATION  OF  THL  W40  COMPLEX 
J.  Smith.  A.  Bintify.  M.  Castilaz.  R  D  Gmirz.  G.  L.  Crash  aits.  and  J.  A.  Hackwti.l 

W  yoming  InfrarcJ  Obstf> jior>.  Departmrni  of  I’hyvies  and  A»ironomv.  University  of  W  yarning 
HtitlttJ /V.fd  August  A  LUxi'PtfJ  /V*4  Oit*»hrr  JO 

ABSTRACT 

Infrared  images  and  photometry  are  presented  for  the  W40  molecular  cloud  and  II  tt  region  complex.  The 
images,  sampled  through  a  5"  aperture  and  observed  Jt  effocme  wavelengths  of  2.3.  3.6.  4.9.  10.0.  and  19.5  pm. 
reveal  onlv  unresolved  sources  clustered  near  the  peak  of  thermal  II  tt  emission  Six  of  the  seven  detected 
sources  coincide  with  faint  but  optically  observable  stjrs.  The  photometry  suggests  that  the  three  brightest 
sources  are  heavily  obscured  early-type  stars  surrounded  by  infrared  emitting  cireumstellar  materials.  In  some 
cases,  much  of  the  compact  infrared  emission  mav  be  produced  by  cireumstellar  condensations. 

The  three  brightest  sources  (IRS  la.  2a.  and  3ui  are  candidates  for  producing  the  level  of  excitation  inferred 
for  both  the  dust  and  H  it  plasma  of  the  W40  complex.  Observed  cireumstellar  dust  has  a  negligible  effect  on 
the  excitation  of  the  W40  region.  Instead,  the  predominant  mode  of  excitation,  whereby  radian!  luminosities  of 
the  exciting  stars  are  absorbed  by  the  \A4t)  dust,  appears  to  involve  only  dilluse  dust  mixed  with  the  molecu¬ 
les.  atoms  and  plasma.  * 

S'lmiecr  headings:  infrared:  sources  —  interstellar:  matter  —  nebulae:  H  it  regions  —  nebulae:  individual 


I.  INTRODUCTION 

W-iOiiar,  H  :t  region  and  molecular  cloud  complex  lying  3.5 
above  the  galjetic  pane  and  approximated  7(>0  pc  !rom  f  arth. 
T  I  it  region  nas  a  diameter  of  -  5  and  a  moueraie  emission 
m — ere  of  -  I  *  10*  pc  cm'*  (Goss  and  Shaver  1970k, 
showing  that  it  has  undergone  the  rapid  phases  of  dynamical 
evolution  normally  associated  with  the  formation  and  expan¬ 
sion  of  a  compact  H  it  region.  Other  components  of  the  W40 
complex  are  a  ! O'  »  20  molecular  cloud  of  moderate  CO 
brightness  temperature  (Zeilik  and  Lada  I9'8).  clouds  of 
obscuring  dust  hat  radiate  bright  far-intrared  and  middle- 
intrarcd  eontinuj  (Oithof  1974:  Price  and  Walker  |9",6i.  and 
a  cTu-ter  of  mtrared  sources  centered  on  the  H  it  region  (Zeilik 
and  Lada  I9'M  Formed  from  these  components,  the  Woo 
morphology  can  be  described  as  an  extensive  molecular 
complex,  the  bneht  core  of  which  lies  in  a  side-ny-side  relation¬ 
ship  with  an  evolved  H  n  region  and  ns  sources  of  excitation. 

In  view  of  its  oistance  from  Earth,  galactic  environment,  and 
observed  morphology,  tf  •’  W40  complex  should  prove  inter¬ 
esting  for  j  variety  ol  detailed  studies  pertinent  to  young  stars 
and  how  thev  interact  wuh  their  interstellar  environments.  For 
this  studs,  wc  chose  to  concentrate  on  a  small.  -3  arcmin* 
region  centered  near  the  position  of  peak  M  n  emission.  Our 
otvective  was  to  produce  a  detailed  infrared  study  of  the  excit¬ 
ing  sources  ot  the  H  ti  region. 

tl.  observations 

All  measurements  were  made  from  the  Wyoming  Infrared 
Observatory  iWIROl.  The  first  step  in  data  acquisition 
involved  taking  64"  <  64”  images  centered  on  the  three  bright 
2.2  urn  sources  discovered  and  named  IRS  I.  2.  and  3  by  Zeilik 
and  '  ada  (I978i  Each  of  the  images  was  obtained  using  the 
mt  is  described  in  Hackwell.  Grasdalen.  and  Gehrz  (19821 
jnd  Uehrz  ft  til.  |  I982l.  For  all  of  the  observed  images,  includ¬ 
ing  those  required  for  the  calibration  stars,  the  sampling  inter¬ 
val  was  I"  with  the  reference  aperture  positioned  60"  north  of 
he  signal  aperture.  The  aperture  diameter  was  —  5'  (FWHM). 


Thermal  infrared  emission  produced  by  the  telescope  and 
sky  must  be  removed  from  the  total  signal  For  the  /.  >  >  6  um 
images,  which  were  observed  with  a  bolometer  detector,  this 
terrestrial  emission  was  corrected  with  a  twrwep  process. 
First,  reference  aperture  signals  were  subtracted  from  the 
image  signals  during  data  acquisition  at  the  W  IRO  Later,  finai 
reduction  involved  removing  from  each  image  anv  etTects  of 
spatial  gradients  in  delected  terrestrial  emission.  These  etTects 
were  interred  from  columns  of  pixels  forming  image  bound¬ 
aries.  along  which  zero  celestial  emission  was  assumej. 

1'nlike  the  /.  S  3.6  inn  images  observed  with  the  bolometer 
AC  coupted  to  its  preampiilter.  the  2.3  am  imaees  were 
observed  with  a  photovoltaic  detector  coupled  directly  in  the 
current  integrating  mode  to  its  preamplifier  and  data  acquisi¬ 
tion  electronics.  Since  reference  aperture  sicnals  were  not  sub¬ 
tracted.  each  pixel  of  integrated  charge  represents  an  absolute 
surface  brightness.  Each  2.3  gm  surface  brightness  was  cor¬ 
rected  for  terrestrial  emission  in  much  the  same  way  described 
above  for  the  /.  >  3  6  urn  images. 

Guided  by  the  preliminary  imaging  results,  we  began  mea¬ 
suring  optical-infrared  energy  distributions  of  the  brightest 
imaged  sources.  Infrared  photometry  was  done  with  reference 
aperture  spactnes  ranging  from  10"  to  60".  For  the  L'BV  pho¬ 
tometry.  the  skv  readmes  were  taken  20"  north  of  each  mea¬ 
sured  star.  Apertures  iFWHM)  were  5"  for  the  4  9  um.  mostly 
5"  but  iKcastonally  10'  for  the  KL,  10"  for  the  JH.  and  12"  for 
the  L  R\  photometry. 

We  calibrated  the  observed  4  9  j<m  flux  densities  and  also 
most  of  the  2.3  and  3  6  dux  densities  with  observations  of 
the  stars  n  LMa.  i  Lyr.  II  Peg,  and  x  Boo  (Gehrz.  Hackwell. 
and  Jones  19741  See  Johnson  et  ill.  ( 19661  for  the  tlux  densities 
of  a  Her  and  v  Cyg.  observations  of  which  calibrated  the  JH 
and  also  the  remaining  K  (2.3  Finland  L  (3.6  /tmi  photometry. 

Positions  of  the  infrared  sources  were  measured  hv  peaking 
up  on  the  individual  sources  and  recording  the  telescope  posi¬ 
tion.  Observations  of  the  IRC  +  00  360,  SAO  142364.  and  SAG 
142368  positions  calibrated  the  W40  positions  (Table  1 1,  each 
measured  with  a  +  1"  uncertainty. 
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tn.  ki.su  rs 

Figure  I  shows  the  five  infrared  images  observed  for  por¬ 
tions  ot  the  V 4<i  complex.  Dashed  lines  give  boundaries  of  the 
observed  images  ranging  in  size  from  5  aremin-  for  2.3  um 
down  to  I  arcmm*  for  4  4  and  14.5  um  See  the  upper  left 
corner  of  Figure  lu  for  representative  point  source  responses 
for  the  photometers  used  to  measure  the  images  E  ach  response 
has  -  2n  arxec*'  for  its  integrated  area  and  F  A\  HM  ^  5'  For 
each  image,  labie  1  gives  the  surface  brightnesses  correspond¬ 
ing  to  the  maximum  contour  levei  of  100  and  to  the  faintest 
level  that  was  chosen  to  have  the  3  rr  statistical  signinearce. 
The  faintest  ievels  arc  given  in  units  of  mJy  arcscc  *.  percent  of 
mjximum  surface  brightness,  and  also  magnitude  per  5  aper¬ 
ture 

The  2.3  *. m  imaec  (F’lg  lui  consists  of  seven  sources  iIRS  la. 
lb.  K.  Id.  2a.  2b.  and  3ai  clustered  near  the  position  of  peak. 
H  n  cmissK'n  labeled  by  "FI  it”  in  Figure  lu  F  or  /.  >  2.3  urn 
if  igs  lu- let.  each  delected  source  hjs  FAVHM  %  5  .  the  size 
•  pected  for  a  point  source  observed  with  our  photometers:  we 
id  no  significant  evidence  fc'r  extended  emission  at  any 
observed  wavelength  The  positions  of  the  infrared  sources  are 
given  in  Table  2.  The  positionsof  IRS.lj.  2a.  and  3a  jgree  with 
the  positions  given  bv  Zciiik  and  Lada  (I4?8i  for  IRS  I.  2. 
and  3. 

A  comparison  of  the  2  3  um  image  with  the  POSS  F.  print 
I  Fig  l/i  shows  that  all  but  one  (IRS  lb)  of  the  seven  infrared 
sources  has  an  optical  counterpart.  Anv  optical  counterpart  to 
IKS  lb  must  have  R  2  20.  the  limiting  magnitude  of  the  E 
print  Figure  lu  and  Tjble  2  give  the  positions  and  names  of 
the  optical  sources  Two  of  the  compact  v'ptical  sources  in 
Tjbic  2.  OS  2c  and  4a.  are  visible  in  Figure  If.  but  we  did  not 
detect  infrjred  counterparts.  Onl>  OS  la.  2a.  and  3a  are  visible 
on  the  blue  POSS  print. 

Of  the  infrared  observations  presented  in  Figure  I.  the 
2.3  // m  image  has  the  greatest  sensitivity  to  main-sequence 
stjrs  obscured  by  moderate  column  densities  of  dust.  The 
10.5  mag  limit  corresponds  to  the  2.3  um  mjg  of  jn  unob- 
seured  B5  ZAMS  star  located  at  the  W40  distance  of  TOO  pc 
Isee  below).  Therefore,  any  main-sequence  stars  too  faint  to 
appear  in  the  2.3  /rm  image  should  contribute  negligibly  to  the 
04  V  level  of  H  n  excitation  inferred  below  for  \V  40. 

The  mean  surface  brightness  of  H  it  continuum  emission 
evtrapoljtcd  from  the  radio  continuum  measurements  of  W40 
(Civiss  and  Shjvcr  1470.  K’Jm  =  4',  Zciiik  and  l.ada  1478. 
beam  =  2T  ts  approximately  20  times  smaller  than  the  limiting 
3  a  surface  brightness  at  2.3  /im  of  •>•2  *  10  ■*  Jy  arcscc  . 
That  limit  and  also  the  small  M"  reference  aperture  spacing 
"•  'Main  why  Figure  lu  shows  no  diffuse  II  ii  emission. 

he  2  2  ;<m  map  made  bv  Zciiik  and  Lada  t  I4'8t  with  Urge 
sizes  of  aperture  diameter  (44  |  and  reference  aperture  spacing 


(240' i  covers  approximately  the  same  region  shown  ir.  Figure 
lu.  Summing  the  2.3  urn  tlux  densities  of  all  seven  infrared 
sources  of  Figure  lu  gives  '.I  Jy.  which  approximate'  closets 
the  10;  1  Jy  (corrected  slightly  from  2.2  to  2.3  pmi  given  bv 
Zetlik  and  Lada  for  their  map  This  means  that  an\  ditlusc 
emission  detected  bv  Zeilik  and  Lada  must  contribute  ie>»  than 
—  3  Jy  to  the  total  2.3  um  flux  density  of  Vk  4a ).  The  tact  mat  the 
contour  map  presented  by  7cilik  and  Lada  appears  to  oc  a  low 
resolution  version  of  Future  lu  provides  turtlier  support  lor  the 
view  tnat  emission  from  compact  inirared  sources  i«  larger 
than  ditlusc  2  um  emission  from  H  it  eas  or  hot  dust  in  Vs -to 
Table  3  gives  the  photometrs  plotted  in  Figure  2  to*  IRS  la. 
lb.  le.  2a.  and  3a  Entries  in  T  able  3  are  in  magnitudes  E  xcept 
where  noted  in  Tabie  3.  I  rr  uncertainties  are  -  U  !  mag 

3  MILL  : 
wa.  ■  Positions 


Names  jiWSh  ivil-sn, 

A  C  .impact  Ss.urcr* 

iks  ia  . .  i«*:s“si-r  .  j  ,  *  xr- 

os  la .  |X  ;s  51  «  <•*  35 

tks  ih .  i«  >  52  6  -7  i  -  a: 

Iks  ic .  b>  4i9  _ 2  r- yd 

os  ic .  lx  49  7  -’I’M 

IRs  id .  lo  ;s  51  4  -2  o*M 

os  id .  is  :x  51  a  - y  o’ 

IRS  2a .  IK  :•»  47  X  -  2  «•*  41 

OS  2a .  IX  2*  47.5  -2  I  *40 

Iks  20 .  JX  2S  4(1  5  -  :  -  45 

OS  2b .  1X2X46  1  -2  "*47 

OS  2c .  IX  2X  490  -2  ti' 21 

IRS  .'a  .  IX  2X  47  H  -2  or.  21 

OS  3a .  IK  :x  47  8  -2  t».  22 

OS  4a  .  IK  28  50  J  -2  10.42 

B.  Large-Scale  Sources 

It  it  .  IK  28  410  07  35 

CO  Peak  t  .  |X  ;x  id  7  -2  O'  5X 

CO-Peak  2  .  |X  2X  39  -  2  isi  25 

ALOL2I77 .  18  2X  47  -  2  0*36 


Non  -  Where  appiopnaic.  an  inirared  source  lIRSi  i\  paired  w.ih  ils 
corresponding  i.pncJl  v.urce  |t)Sl  Inirared  source  p.silions  am-  mea¬ 
sured  al  2  3.  3 6,  or  l.i  um  and  w.ih  5  apertures  Optical  source  p'silions 
were  measured  on  a  print  ot  ihc  Palomat  I  Ibsen  alors  Sk  v  Survec  t  -arrrev 
a  pinion  ol  which  rs  reproduced  in  I  ig  1/  Positions  lib.  telerring  In 
large-scale  comp.nrnis  of  Ihr  SX40  complex  were  measured  wilt  -3 
apenures  Pubrohed  I  rr  uncertainties  are  -  20"  lor  the  p.stlion  o'  peak 
It  it  emission  lOoss  and  Shu.er  |9**ii»  rfnd  -  2  for  Al  (>l  21"*  iPrts'  jnd 
Walker  I976r  Pool  ions  of  peak  *(  OlZ  ■  I  *  OS  emission  were  taken 
from  /riltk  and  Lada  I97g. 
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TABLE  3 


WAI  Photomi  Tin 


3  Source 

IRS  U 

IRS  16 

IRS  lc 

IRS  2a 

IRS  3a 

t  -« . 

Oh 

1  1 

OV 

« - »  . 

2  "v 

2  7 

24 

1  . 

ISO 

16  1 

15  6 

j  a . 

1  3 

1  6 

09 

II  -  K  . 

1  0 

09 

0  5 

K-L . 

1  3 

20 

It 

i '  i  o .: 

04 

1  . 

4  J 

5  7 

4h 

-  \ 

L-  \t . 

Oh 

0* 

O' 

t  ii 

O' 

i  •  |s  '  i.m J . 

1 ' 

24 

22 

3* 

2."  -  0.2 

L-S  . 

1  s 

^  2 

:  4 

4  5 

3.: 

L  -  [  1 1  4  i.m]  . 

2.2 

3  1 

2.6 

JO 

3  4 

L  -  I  12  n  um  j  . 

2.3 

:  v  -  o.: 

53 

4  o  r  0.2 

L  —  I  ?9  5  um]  . 

j.: 

4  “  -  04 

6' 

5.3  *04 

1  <m|  . 

7.5  ♦  0  2 

The  shapes  of  the 

infrared  energv  distributions 

given  in 

Future  2  demonstrate 

that  the  three  sources  seen  in  th 

e  l<>  i<m 

image*  (IRS  la.  lb.  and  2.n  have  eircumsteiiar  emission  Fiux 
densities  of  IKS  2j  increase  from  lit  to  23  i/ni  in  j  tjshion 
expected  of  warm  eircumsteiiar  dust  Infrared  sources  la  anj 
Ih  shou  broad.  roughly  constant  distributions  of  intrared  Dux 
dcnsits  possibix  from  eircumsteiiar  material  having  a  wide 
temperature  ranee  Notice  tnat  pnotospnenc  radiation  from  a 
main -sea uence  star  would  haxe  a  monotomcaliy  decreasing 
flux  density  as  the  infrared  wavelength  increases. 

Figure  2  presents  exidence  for  eircumsteiiar  emission  from 
two  additional  sources  <  1 R S  lc.  3ai  too  fan  t  to  be  imaged  at 

i  nm  At  /-  2  I*1  uni.  IRS  3a  resembles  IRS  2a.  whereas  IRS 
.c  has  a  broad  infrared  energy  distrinution  much  lice  that 
shown  for  IRS  la  Tne  small  amount  of  information  given  in 
Figure  2b  for  IRS  Id  and  2b  was  inferred  from  the  2.3  and 
3  b  urn  images  of  Figure  1. 

F  rom  3  b  to  12b  cm.  the  cncrgv  distribution  shapes  iFig  26) 
of  IRS  In  and  lc  resemble  the  approximates  Slat  energx  dis¬ 
tribution  of  their  nearest  bright  neighbor.  IRS  la  A  note- 
worths  difference  among  the  three  energy  distributions  is  the 
abrupt  decrease  occurring  between  3  6  and  2.3  pm  for  IRS  lb. 
the  only  source  of  the  three  that  has  no  xisible  companion  on 
the  ROSS  E  print.  Possible  explanations  of  these  two  differ¬ 
ences  in  IRS  lb  arc  that  it  has  a  larger  extinction  or  an  intrinsi¬ 
cally  fainter  star 

7eilik  and  Lada  <I97S|  give  some  energy  distribution  infor¬ 
mation  for  their  sources.  IRS  I,  2.  and  3  Their  data,  obserxed 
with  three  different  aperture  sizes  and  a  reference  aperture 
spacing  of  100".  can  be  compared  wuh  our  detailed  energy 
distributions  of  IRS  la.  2a,  and  3a.  each  observed  with  a  5" 
aperture  and  a  60  reference  aperture  spacing  Zeilik  and  Lada 
observed  IRS  I  with  an  IK  aperture  At  K  and  L.  their  photo¬ 
metry  of  IRS  I  and  our  photometry  of  IRS  la  agree  reasonably 
well.  Beyond  3  6  pm.  the  photometry  of  IRS  I  and  IRS  la 
differ  the  energy  distribution  of  IRS  I  increases  strongly  from 
5  to  20  pm  in  a  fashion  expected  of  wjrm  dust,  whereas  IRS  la 
shows  a  broad  and  nearly  llat  energy  distribution  The  energy 
distributions  given  by  Zeilik  and  Lada  for  IRS  2  and  3  cover 
the  small  wavelength  ranges  of  165-5.0  and  1.25-50  pm. 
respectively.  Observed  through  a  16"  aperture,  their  photo- 
ctry  ol  IRS  2  agrees  reasonably  well  with  our  photometry  of 
<S  2a  In  contrast,  their  photometry  of  IRS  3.  observed 
through  a  large  46”  aperture,  differs  from  our  photometry  of 


IRS  3a :  IRS  3  is  fainter  at  J  and  H  but  brighter  at  the  larger  K 
and  L  wavelengths  In  summary,  the  photometric  observations 
reported  by  Zeilik  and  Lada  are  brighter  than  or  equal  to  our 
values  at  K  and  longer  wavelengths:  at  the  shorter  J  and  H 
wavelengths,  their  datj  ate  euual  to  or  fainter  than  ours 
The  way  in  which  the  VvIRO  and  Zeilik  and  Lada  energy- 
distributions  depend  on  aperture  size  and  reference  aperture 
spacing  suggests  that  W  -40  has  extended  distributions  of  infra¬ 
red  emission.  Some  good  evidence  for  extended  emission  from 
interstellar  dust  is  the  large  ratio  of  IRS  I  IRS  la  Dux  densities 


10  0 


X  l  p  m  ) 
3  C 


0  3 


o)  4  i«s  la 
o  IRS  2  o 
•  iRS  3  a 


w 

\t 


LOG  v 

flci  ?  -  \  nerpv  distributions  of  thf  soen  infrared  source*  observed  in 
W40  Logarithm  of  flu*  dcnsits  if,  in  J>i  is  pinned  npjinvt  log.icthm  of 
radiation  frevjuenc*  Solid  line*  show  smooth  trend*  in  the  data  plotted  for 
ejLh  source  Set  1  ihle  3  lor  the  broaU-banJ  photometr>  plotted  for  IRS  la, 

lb.  lc.  -a.  and  >j  f  lu*  densities  of  IRS  Id  and  2b  Here  inferred  !r««m  the 
infrjred  images  shown  in  F  igx  la  and  lb  tui  Optical  through  infrared  energy 
disinbutions  of  IRS  la.  2j.  jnd  ^a  tM  Infrared  enerex  distributions  ol  IRS  lb. 

lc.  Id.  and  2b  The  upper  limit  plotted  at  3  6  Mm  for  IKS  2b  is  3  e. 
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IR  SOLRCLS  AND  LXCITATION  OF  W'40  COMPLEX 
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X(^m) 


LOS  9 

f'li 1  — Infrared  energy  tJi>inf*uti“ns  of  the  VlJj)  ^orrmex  Lofjnthm  of 
tr.e  pr*HJw^.i  of  frcuutnc*  timex  oen-*it>  t.f in  \fc  m  *»  ploitea  jiainst 
lotfjrnr.rr.  i  the  treuucno  l  uendinii  Irom  1.3  Jo  Is*  5  -m.  the  *on«l  '.ne 
(he  rie-in*.  .Jjt  'pcctrjt  di>trir>utn>n  «*l  the  >um  kwh  .  iti  ini  ol  the  ilu*  .ien\itie* 
mea%oreo  :»»r  tnense  nrter.icM  ifiirjrevJ  tourers  ami  w:*en  inuowluailv  in  fjNe 
3  l  or.nr.'ea  b*  the  ca>neu  hnc  t.iied  -ircies  xn.**  the  •*  U.  -and  20  *(n 
rrcjsurer'sm  made  Prve  and  t;ker  «S^>i  v^-tn  a  3  5  «  10  5  jrxrture 
and  tr.e  ar.d  155  um  measurement*  made  nv  Oi!otiW*4i  *nh  a  ;0  aper¬ 
ture  The  m.»notonicjtU  tm-rca%ifliz  dotted  line  t>  an  une\t»neui*fteu  eiiraptw 
Uiton  ol  -he  i  4-  5.0  OH*  H  it  emivuon*  measured  b>  Vtennoif  ei  M  « W“0i 
*ur*  an  1 1  aperture. 


amounting  to  j  factor  of  -10  in  the  10-20  i-m  intenal  In 
addition.  pnotometnc  errors  mas  contribute  to  the  laree  differ¬ 
ence  between  observations  of  IRS  1  and  la.  The  differences 
observed  between  IRS  3  and  IRS  ?a  ma>  reouire  extended 
emissions  of  Iwo  types:  emission  from  hot  interstellar  dust  mas 
explain  tne  differences  observed  at  K  and  L.  whereas  photo- 
spheric  emission  produced  by  the  widely  distributed  infrared 
cluster  >f  W4<)  iSmith  19V4t  max  explain  the  differences 
observer  at  the  shorter  J  and  H  wavelengths. 

The  presentation  of  Figure  3  reveals  some  additional  infor¬ 
mation  concerning  the  origin  of  infrared  emissions  observed  in 
Wit).  Total  flux  densities  were  inferred  irom  the  large  beam 
measurements  of  Price  and  Walker  11976)  and  Olthof  (I9"4i 
made  at  middlc-inlrared  and  far-mlrarcd  wavelengths,  respec¬ 
tively  The  contribution  ol  compact  infrjred  sources  to  each 
total  ffux  density  is  approximated  by  the  sum  ol  the  tlux  den¬ 
sities  measured  for  IRS  la.  lb.  !c.  2a.  and  3a  (Table  3);  using 
the  low -resolution  observations  of  IRS  I.  2.  and  J  made  by 
Zcilik  and  Lada  ( 197X1  would  make  no  qualitative  difference  in 
our  discussion.  For  observed  wavelengths  of  the  10-20  um 
interval,  the  combined  emission  of  the  live  infrared  sources 
accounts  for  <  10"  .  of  the  total  emission.  This  small  contribu¬ 
tion  of  the  infrared  sources  and  also  the  shape  of  the  W'40 
energy  distribution  suggest  that  diffuse  clouds  ot  dust  having  a 
wide  temperature  range  produce  almost  all  the  middle-infrared 
and  far-infrared  emissions  observed  for  W40.  About  9b'’..  of 
the  total  W'40  emission  is  produced  beyond  20  /<m,  suggesting 


that  the  cool  component  of  the  diffuse  dust  is  the  predominant 
source  of  reradution  in  W'40.  Much  of  the  4  remainder 
emitted  short  of  20  iim  jppears  to  be  produced  by  a  warm 
component  of  the  diffuse  dust  possibly  closely  associated  with 
the  M  it  plasma.  Near  4  um.  the  total  emission  must  be  produc¬ 
ed  by  a  complex  combination  of  sources  includific  the  infrared 
sources  and  the  H  n  gas.  each  of  which  contributes  -25"..  of 
the  total  amount. 

Before  the  levels  of  excitation  can  be  estimated  for  the  diffuse 
clouds  of  dust  and  If  it  gas  in  W'40.  their  distance  must  he 
specified.  The  distance  adopted  for  our  work  is  "00  pc  iGoss 
and  Shaver  1970).  This  value  was  inferred  from  recombination 
line  observations  of  the  W40  H  it  region  and  the  kinematic 
model  of  the  Galaxy.  Another  datum  that  tends  to  eonlirm  this 
relatively  small  distance  of  W'4t)  is  as  galactic  latitude  of  3-5. 
This  latitude  requires  that  W'40  be  no  more  distant  than  -  MX) 
pc  if  it  is  to  reside  within  the  -  kx)  pc  thick  plane  normally 
populated  by  the  Galaxy  's  H  it  regions. 

The  level  of  dust  cloud  excitation  may  be  inferred  from  the 
infrared  observations  given  for  the  W 40  complex  in  Figure.  3. 
Summing  the  fluxes  implied  by  the  dashed  curve  gives 
5.x  10  "*  W"  for  the  4-21X)  um  flux  of  the  W'40  complex. 
That  flux,  or  its  equivalent  luminosity  of  "5  *  lO*  L  .  could 
ortemate  from  dust  heated  by  a  single  09  V  star.  This  same 
09  \  level  of  excitation  can  be  calculated  for  the  M  it  eas  from 
radio  observations.  For  the  H  u  region  of  W'40.  Altenhoff  ct  til. 
(|9iWi  give  35  Jv  for  the  radio  continuum  emission  and  Pank- 
onin.  Thomassohn.  and  Barsuhn  1 1977)  give  He  '  H‘  =0.09 
for  the  lisisx  ratio  of  recombination  line  fluxes.  Both  the 
Lvnjn-continuum  luminosity  of  -15  *  Hr*’  photon  s  ’ 
inferred  from  the  continuum  emission  and  the  approximate 
excitation  temperature  implied  by  the  recombination  ratio 
could  be  produced  by  a  single  09  \  star. 

IV.  DIXLSSION 

ui  Fneruy  Distributions  ot  the  ('ompiicf  .Votaries 

The  evaluation  of  jeeurate  momxihromatie  extinctions  is  a 
crucial  part  of  any  attempt  to  decompose  jn  observed  energy 
distribution  into  ctrcumstcll.tr  and  pnoiosphene  components. 
Table  4  cives  results  ot  the  sequence  of  steps  leading  to  the 
mom>ehrom.itic  extinctions  ot  W4"-IRS  :a.  fust,  we  lounj 
that  the  IRS  3a  photometry  of  the  •)  55-49  „m  interval  could 
be  explained  by  the  photospherte  emission  of  a  Bl  V  star  lying 
at  the  W  4t)  distance  of  "IX)  pc  and  obscured  by  a  standard  tv  pc 
of  interstellar  extinction  (Savage  and  Mathis  |9"9l.  As  a  result, 
the  Fit'  -  /. I  excess  colors  were  compuls'd  as  the  observed 
colors  (the  F  -  [/.]i  minus  the  intrinsic  colors  of  a  Bl  V  star 
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(Johnson  1 9661.  Using  the  five  excess  colors  of  (he  0.55-4  9  pm 
__  interval  and  standard  interstellar  values  of  A,  hit'  —  /.I.  we 
pmputed  100  mac  for  the  mean  visual  extinction  of  IRS  3a. 

*  t  he  last  step  of  the  extinction  analysis  was  to  use  the  derived 
10.0  mag  of  visual  extinction  to  compute  each  of  the  values  of 
monochromatic  extinction  (.4,1.  for  /.  £  0.55  pm.  each  value  of 
.4,  was  computed  as  the  product  of  the  10.0  mag  and  the 
appropriate  standard  value  of .(,  A,  For  the  shorter  H  and  L' 
wavelengths,  the  monochromatic  extinctions  were  computed 
directlv  from  the  excess  colors  as  I0.H  Li/l-Ci  and  10.0 
■*  HI  -  •  I  Values  of  .4,  f.lH-  IT  and  4,  El f-  I  t  implied 
by  this  anjlvsis  are  3  h  and  2.2.  -respectively  Diflcnng  little 
from  the  standard  xjluescf  3  I  and  1.9.  those  extinction  ratn>s 
emphasi/e  the  predominantly  standard  nature  of  the  mono¬ 
chromatic  extinctions  aliectine  the  observed  energy  distnou- 
tion  of  IKS  3a 

Figure  4  shows  the  energx  distribution  of  \V40-IRS  3a 
decomposed  into  phoiospheric  and  circumsieliar  components. 
The  pholosphenc  energy  distribution  is  normalized  to  the  I 
flux  density  of  IRS  3a  and  has  the  shape  determined  bv  the 
intrinsic  colors  of  a  Bl  V  star  modilied  bv  the  monochromatic 
extinctions  of  TaPie  4  The  circunistellar  energv  distribution  of 
IRS  3a.  derived  bv  subtracting  the  photospneric  (lux  densities 
from  the  observed  values,  is  indicative  of  warm  circumsieliar 
Uusi  at  -s  2"’i|  K  The  circumsieliar  luminosity  of  the  5—2*  1  urn 
interval  is  4  L  or  -  U  OF  ..  ot  the  total  luminosity  of  a  B I  V 
star 

L  nlike  the  colors  of  IRS  3a.  the  colors  of  IRS  la  and  2a  do 
not  show  convincing  evidence  for  predominantly  photospneric 
emission  extending  from  optical  to  ncar-inlrared  wavelengths.  - 
Therelore.  we  must  seek  an  alternative  to  tne  method  of  inter- 
ng  a  mean  value  of  visual  exiinction  from  the  observed 
opticai-imrared  co  srs.  A  reasonable  alternative  is  to  use  the 
IRs  3a  ratios  of  4,  Fifl— f  i  and  .4,  f.iL  —  IT  together  with 
values  of  the  excess  colors.  ElB-  IT  and  Ft l  -  1  I  For  each 
star,  the  excess  colors  were  computed  from  the  observed  colors 
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ft«*  *  - haiirnisiMiion  of  the  W40IKS  3a  energv  distribution  inio  cir- 
cumsirlljr  am)  pholosphenc  components  Observed  data  are  plotted  as  filled 
voiles  Data  observed  lor  ihe  ll.sX-JVpm  interval  are  lined  to  the  phoio- 
P  eric  emission  of  a  bl  V  phoiosphere  tJa\h  J  linn  obscurrd  bs  iniersiellar 
usi  producing  |.  -  III  mag  Shown  connecled  bv  the  wild  line.  Ihe  urcum- 
Vlellar  lluv  demines  are  the  observed  minus  ihe  pholosphenc  flu»  densities. 
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F io  5  — Fnerpv  distributions  r»f  Wan  sources  corrected  for  rp'ersicllar 
extinction  Lash  dolled  line  shows  me  share  ol  a  HI  V  enerpv  di-iribulion 
nvirman/cd  lo  me  corrected  I  ftuvdenvilv.  4  displav  consiam  i(  lajcrdtoihe 
loeariihTn  of  lluv  uensns  if.  in  Jsi  is  pioued  aeainsi  logarithm  .  f  rad. anon 
(reuuencs  Values  of  C  are  +3.  -4  3.  and  -4  lor  IKS  la.  da.  and  la.  respec¬ 
tively 

(Tabic  3) and  Bl  V  as  a  nonen tical  approximation  to  (he  spec¬ 
tral  type.  Values  of  mean  visual  extinction  inferred  by  this 
alternative  method  arc  9.2  and  11.3  mag  for  IRS  la  and  2a. 
respectively.  Monochromatic  extinctions  are  inferred  from 
each  of  these  values  of  visual  extinction  in  the  same  way 
described  above  for  IRS  3a. 

Figure  5  shows  the  energy  distributions  of  IRS  la.  2a.  and 
also  3a  corrected  for  the  interstellar  extinctions  computed 
above  Also  plotted  in  Figure  5  arc  the  energy  distributions  of  a 
Bl  V  photosphere  normalized  to  the  corrected  visual  flux 
density  of  each  star.  Although  circumsieliar  emission  becomes 
important  for  IRS  3a  only  at  wavelengths  longer  than  4  9  pm. 
it  dominates  the  energy  distributions  of  IRS  fa  and  2a  at  much 
shorter  wavelengths. 

It  proves  informative  to  compare  the  W40  sources  with 
other  objects  w  hich  have  been  studied  in  some  detail,  F  igure  6 a 
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Ft*#  p  -Encrr.  distribution*  *>t  cornnanvn  >urs  uv*U  t*»  Muds  the  oPvrvcd  '«U*i  source*  F»>r  each  star  v'lihc  two  panels  j  d;spuv  const  am  ii’i  joded  to  the 
loiunihr?  of  riu*  census  »F  .  n  J v » t-»  pioticd  .ersus  the  toijrttP.m  oi  on^erxed  radutum  regueno  t.ii  A  -eguence  t*l  hot.  cmiwion-une  Mjr\  ordered  .»aord:na  to 
their  Id*  «)  >5  am  ratio  oi  flu*  densities  3imuar  u>  \A4iwlRs  la.  lb.  and  iv.f  ,  *  wonoart  lor  *  uPxianiijI  pavilion*  «»l'  the  infrared  enerss  distributions  v>i,d  ojfin 
shot*  mereis  »rends  in  me  rtu«  uenut  et  Values  oi  C*  are  ♦  3.  0,  -  *.  and  -~n  for  T  On  iB*  A  jr*i.  V  'ho  Dn'B>  A2e».  R  Mont  -e  ♦•»nciU  ind  V  P-i*  c’.ani  *  \5<l 
respective  v  References  lor  the  Plotted  rt us  densities  are  VA  alker  i  Pension  *  I'**’  t».  and  Ntrith  i  W**»  lor  f  On.  Lee « tM'Oi  Loss  er  j..  « 4  or.en  i  and 
Alien  for  v  •sit On  and  R  Mon.  and  H  irrpnresx  Mernil.  and  Black  t  tor  VMM  v am  References  for  the  emisuon-une  Npcvira  are  J»  nnson  i  :•**<•  \»r  T 

On  HerNa  il^^O.  ; *#0  for  V  *M>  t>n.  and  R  Mon.  Humphress  er  at  lor  VM'  i  ,*m  »ni  |  xtm«.tion  seuuence  ot  compartson  OB  'tars  re:at.*d  m  Wai'-JRS  '*x 
Svmrmi>  m  the  oosersed  t!uv  densities.  Photo>rncnc  tius  densities  iJu.wirj  /inr»i  *ere  computed  tr»»m  intrinsic  stcilar  coiors  miniiPed  bv  ;n*erstei!ar  esttn^iioit 
Connected  bv  a  solid  line  for  each  star,  ^ircumsiei.ar  riux  densities  are  the  orservet]  minus  me  pnotosphenc  rtu\  densities  V  alues  oft  4  .  O  are  »'*.  ft.  1 1  I.  -  ' <* L  *2  2. 
-  -4ti».  ard  i*s  maa.  -POi  tor  tne  intrinsic  energv  distnbuiion  ot  a  B0  V  star  nop.  V  «  o  maat.  t‘*A  On  iO*F5  V  l  NL  On  tB0.5  Vi.  and  NtlC*  -0^4  No  l  t  bottom. 
BO  5  Vi.  respectively. 


shows  the  energy  distnbutions  of  four  comparison  stars.  T  On. 
V3N0  Ori.  R  Mon.  and  VM5  Cvgnt.  each  of  which  is  a  hot 
emivMon-line  star  associated  with  nebulosity.  Like  W40-IRS 
la.  lb.  and  lc,  the  emission-line  stars  tend  to  have  F,  % 
constant  tor  a  substantial  portion  of  their  2-20  um  energy 
distributions.  Thus.  IRS  la.  lb,  and  lc  may  be  hot  emission- 
line  siars.  f  igure  bb  shows  the  energy  disiributions  of  three 
additional  comparison  stars,  0:  -\  Ori.  NL'  On.  and  NGC  2024 
No.  I.  each  ot  which  is  a  relatively  normal  OB  star  with  a 
reradiating  circumstellar  dust  cloud.  Like  W4IVIRS  3a.  each  of 
these  OH  stars  has  a  circumstellar  energy  distribution  that 
increases  strongly  from  5  to  20  urn.  Thus,  the  case  is  strength¬ 
ened  for  IRS  3a  being  a  main-sequence  OB  star  surrounded  by 
a  eloud  of  reradiating  dust  grains.  Figure  3  shows  that 
W40-IRS  2a  has  characteristics  of  both  the  emission-line  and 
main-sequence  OB  stars:  for  the  1-5  *im  interval,  the  corrected 
flux  densities  of  IRS  2a  are  nearly  constant,  whereas,  for  wave- 


leneths  longer  than  5  urn.  the  energy  distribution  of  IRS  2a 
increases  strongly  like  those  ol  the  main-sequence  OB  stars. 

The  summary  presented  in  Table  5  includes  descriptions  of 
W40-IRS  3a  and  the  three  comparison  stars.  \  Ori.  NL'  On. 
and  NGC  2024  No  1.  The  above  jnalvsis  of  IRS  3a  gives  the 
guidelines  used  to  derive  monochromjue  extinctions  for  each 
of  the  comparison  stars.  For  eath  of  IRS  3a.  0‘  \  Ori.  NL'  Ori. 
and  NGC  2024  No.  I.  the  value  of  Lr  L,t  shows  that  circum- 
stcllar  emission  I /..  I  is  a  substantial  portion  of  the  total  infrared 
emission  tLwi.  However,  reiaiive  to  the  photospheric  emission 
If.. I  estimated  for  each  star,  the  value  off.,  is  small.  Thus,  very 
little  of  ejeh  star's  luminosity  is  reradiaied  bv  circumstellar 
dust  at  5-20  fim. 

Because  the  photospheric  components  are  very  difficult  to 
identify  for  IRS  la  and  2a,  the  summaries  presented  in  Table  5 
for  these  two  stars  are  relatively  uncertain.  Nevertheless,  these 
summaries  should  be  useful  for  purposes  of  comparison.  Each 
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of  the  three  sources.  IRS  la.  2<t.  and  5a.  consists  of  a  heasily 
upscdied  OB  star  associates]  with  a  compact  circumstellaf 
source  This  interpretation  differs  in  two  basic  ways  from  the  _ 
stew  developed  ns  Zeiiik  and  Lada  iI9“m  for  their  low- 
—solution  opsersations  of  W40.  T'nes  dcscnned  their  th-ee 
.tghtesi  sources  not  as  being  compact  but  as  resolsed.  -  50" 
diameter  regions  of  diffuse  emission  Second.  the>  argued 
that  extinction  h>  dust  should  not  atfect  substantial!)  the  \V4(i 
energy  distributions.  whereas  our  analysis  implies  suost.init.il 
saiues  of  -  10  mac  for  the  visual  extinctions. 

hi  .Morphology  anJ  Lnerueius 

Basic  elements  of  the  W40  morphology  tFic  7|  include  <  I  f 
the  brightest  an.!  therefore  probably  densest  portion  of  the 
extensive  molecular  complex  observed  by  Zcilik  and  Lada 
I  |97X|.  i2i  the  large.  -5  diameter  H  it  region  i  Cross  and  Shaver 
1971.';  Mter.hotl  rt  a/  1970i  King  in  a  sidc-bv-side  and  possibly 
contiguous  relationship  with  the  bright  molecular  core,  and  <51 
the  seven  detected  infrared  sources  distributed  within  a  !!5 
diameter  region  centered  approximately  on  the  position  of 
peak  II  II  emission. 

The  close  association  of  the  infrared  sources  with  the  H  it 
region,  the  absolute  visual  magnitudes  (4/,|  summarized  tn 
Table  5.  and  the  level  of  excitation  derived  for  both  the  dust 
and  H  n  gas  <i  111)  suggest  that  some  combination  of  the 
brightest  detected  sources  iIRS  la.  2a.  5ai  has  enough  lumin¬ 
osity  to  ionize  the  II  it  gas  and  heat  much  of  the  dust  of  the 
W40  complex.  IRS  2a  appears  to  be  the  most  luminous  star; 
we  suggest  that  it  is  the  most  important  single  source  of  excita¬ 
tion. 

Estimates  of  extinction  optical  depths  constrain  the  possible 
modes  whereby  the  obsersed  sources  heat  the  dust  distributed 
throughout  the  W40  complex.  To  compute  extinction  optical 
.  "'cpths  for  the  dust,  we  assume  that  the  dust  extinctions  obey  a 
andard  interstellar  curse  iSasage  and  Mathis  I979|  and  that 
optical  reddening  is  related  to  the  nucleon  column  density 


according  to  the  standard  interstellar  relation  rven  by  Bohim. 
Savage,  and  Drake  t!978i.  These  two  assumptions  lead  to 
s  I  •  10  VM  mag.  the  desired  relation  fviween  the 
ultraviolet  extinction  I 1  and  the  nucieoti  column  density 
t.\Mi.  L  sing  this  relation,  we  estimate  ultraviolet  extinctions  o' 
<5  and  -  5o  mag  for  tne  columns  of  dust  associated  wuh  the 
H  u  and  molecular  gases,  respectively.  For  the  H  t:  region.  \H 


A  R  A  .  (  oremir.) 

T Hi.  ?  —  fclermrn t*  of  the  W40  rorrrhoIoRy  uherr  tmission 

from  the  CO  »nuictulc  hds  been  itctcvie-J  mihc  14  *  1 4  rzfion  cmcrcd  the 
figure  The  v*hd  hnc  envloves  ihe  bright  molecular  ct»re  ior  ^hu.h  the  peit 
bnifhtness  temper  jturo  »•(  CO  \J  ■  I  *  Oi  emisvion  are  »r  ihe  15  31  k  rjnee. 
The  pmiuon*  t>i  pe^k  CO  rmiiwi'n  are  plotted  a>  filled  tramrle^  C  entered  at 
the  tn»*.  plotted  -*2  5  northwest  of  the  poMiion  of  bri^lnw  CO  cmixMon.  the 
5  dumeier  H  u  rranm  /inn  lies  vide  b>  vide  uitft  v*e  mx'levular  ct^re. 

T died  ctrclev  »hoiv  posmonv  of  the  seven  detected  inlrarea  vmrvcv  all  ot  ^hich 
he  nr ar  the  p**viti»»n  of  peak  H  it  emiwi%»n  Povitums  in  arc  minutes  are  reck* 
oned  relative  to  the  IRS  2a  position  plotted  as  <0.  Ot.  At  the  adopted  distance  of 
7du pi.  the  scale  is  I  *l)2opc 
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was  equated  lo  the  »l  *  I011  on  ' 1  electron  column  density 
(Shaver  and  Goss  19701;  .4„,  is  an  upper  limit  because  dust  that 
would  ordinarily  be  associated  with  the  -  I  <  10* *  cm  2 
column  density  of  nucleons  may  be  modified  or  expelled  by  the 
harsh  ultraviolet  radiation  field  of  the  H  u  region.  For  the 
moleeular  complex,  wherein  \M  =*  2.Y,,.  applies,  we  computed 
a  peak  value  of  \M  from  \Ml  =  5  *  itJ'.V,,  (Dickman  I97xi 
and  \(J  =>  5  <  10"’  cm  ' 2.  the  1 JCO  column  density  given  by 
Zeilik  and  Lada  (I97SI  for  the  position  of  Peak  I  (Table  !|. 
Using  the  relative  distribution  of  ,:COiJ  =  I  — » Oy  emission 
(Zeilik  and  Lada  19781  to  scale  the  peak  column  density  gives 
AN  t  J  x  10"  cm  1  for  molecules  distributed  along  the  lines 
of  sight  intercepting  the  embedded  infrared  sources. 

Now  that  we  have  proposed  the  maior  sources  of  excitation 
and  estimated  some  extinction  optical  depths  for  the  volume  of 
Wat),  we  can  proceed  to  discuss  some  general  aspects  of  the 
dust  cloud  energetics.  A  small  portion  of  the  total  amount  of 
ditfuse  dust,  surrounding  the  cluster  and  having  ,4„  %  t.  could 
absorb  and  then  reradiate  muen  of  the  cluster's  luminosity. 
Depending  on  whether  dust  is  depleted  in  the  If  ii  region  and 
also  on  the  accuracy  of  the  above  extinction  calculation,  much 
of  this  dust  with  .4..  s:  I  may  reside  within  an  interlace  • 
between  the  H  n  boundary  and  molecular  cloud.  This  leaves  a 
large  amount  of  cooler  dust,  mixed  with  the  molecules  and 
isolated  by  signiticant  ultraviolet  extinction  from  the  energetic 
sources  of  the  H  u  region,  to  be  heated  by  radiation  ;rom  j 
complex  commnation  of  sources  including  the  onscurgd  excit¬ 
ing  stars,  additional  memoers  of  the  widely  distributed  W40 
star  cluster  (Smith  19X4i.  and  the  warm  dust  grains  near  the 
H  it  region  Each  of  the  two  dust  components  will  contribute 
to  the  total  reradiation  spectrum  of  the  W4()  complex,  but  the 
small  amount  of  diffuse  dust  irradiated  directly  by  the  exciting 
stars  wiil  contribute  most  to  the  far-infrared  rcradiation. 

A  smail  amount  of  ditfuse  dust  may  he  within  the  circumstel- 
lar  region  of  each  of  the  infrared  stars.  Although  this  hvpotheti- 
cal  component  would  have  a  negligible  impact  on  the 
energetics  ot  W  0i).  we  can  show  that  its  presence  eouid  account 
for  a  agnihcant  portion  of  the  ctrcumsteilar  emissions 
observed  for  IRS  2a  and  3a.  We  ur  estimate  the  m/o  of  a  cloud 
oi  ditfuse  dust  for  which  /-  L,  <  t'.i'I.  a  limit  to  the  levels  given 
fv'r  the  7-20  urn  reradiations  of  IRS  2a  and  3a  i  TaPie  >|.  Along 
ns  diameter,  the  cloud  would  produce  4,„  a:  1L.  L.  for  the 
extinction  at  ultraviolet  wjveienctbs  0r  <  I  \.  of  the  ^3  mag 
estimated  for  the  assumed  unmodified  dust  mixed  with  the 
piasma.  If  5 .  roughly  the  observed  width  of  the  H  it  region, 
approximates  the  distance  to  which  the  3  mag  refers,  then  the 
ditluse  cloud  would  be  <2"  in  diameter.  In  addition  to  this 
compact  st/e,  we  would  expect  the  ctrcumsteilar  component  of 
the  ditfu.se  dust  to  have  a  10  20  i.m  color  temperature  larger 
than  the  2Di>  K  interred  from  ihe  AFGL  measurements 
made  with  a  large  beam  iKtg.  3i.  Both  expected  properties  of 
the  ditluse  ctrcumsteilar  emission  are  consistent  with  the  obser¬ 
vations:  the  images  show  that  IRS  2a  is  compact  at  both  it) 
and  20  jim  and  the  sets  of  photometry  give  ^250  K  for  the 
10  20  um  color  temperatures  of  IRS  2a  and  3a. 


The  differences  in  shape  observed  for  the  infrared  energy 
distributions  of  IRS  I  j.  2a.  jnd  3a  suggest  some  intrinsic  differ¬ 
ences  among  (heir  ctrcumsteilar  regions.  Harvey,  Tlironson. 
and  Gatlcy  t!979i  have  discussed  how  the  shape  of  a  circum- 
stellar  energy  distribution  may  depend  on  the  density  protile  of 
the  ctrcumsteilar  material.  One  of  their  models  suegests  that 
strongly  increasing  energy  distributions  of  the  type  observed  at 
/.  >  10  /im  tor  IRS  2a  and  3a  arc  produced  by  circumstelljr 
dust  grains  having  j  uniform  density.  Within  the  context  of  this 
model,  our  above  discussion  implies  that  the  uniformly  distrib¬ 
uted  dust  surrounding  IRS  2a  and  3a  may  be  a  small  portion  of 
the  ditfuse  interstellar  medium  of  WUij.  \  second  model  con¬ 
sidered  by  Harvey.  Thronson.  and  Galley  (1979)  suggests  that 
broad  energy  distributions  of  the  type  observed  for  the 
emission-line  stars  V3S0  On  and  R  Mon  (also  Fig.  6ul  are 
produced  bv  ctrcumsteilar  condensations  of  dust.  \V  4t  I- 1 R S  |j. 
and  possibly  IRS  lb  and  Ic.  have  broad  energy  distributions 
resembling  those  of  the  emission-line  stars  and  therefore  are 
candidates  for  ctrcumsteilar  condensations. 

Being  clustered  tightly  within  a  small.  ^0.2  pc  diameter 
region  of  the  W40  complex,  the  seven  W4()  sources  may  have 
approximately  the  same  epoch  of  formation.  If  so.  these  obiects 
would  have  rougnlv  the  same  age  yet  they  have  a  diverse  set  of 
properties.  For  example.  IRS  la  has  the  F  =:  constant  tvpe  of 
ctrcumsteilar  spectrum,  whereas  IRS  3a  has  ctrcumsteilar  emis¬ 
sion  thjt  increases  strongly  beyond  -  7  urn.  Perhaps  this 
diversity  ot  form  results  from  the  way  that  a  compact  omect 
evoives  within  its  H  u  region  and  molecular  cloud  complex.  In 
the  ease  of  W4().  the  evolution  appears  to  depend  in  a  eompiex 
wav  on  the  individual  evolutions  not  only  of  the  err.rcdued  star 
and  jnv  ctrcumsteilar  condensation  it  mav  have,  but  also  the 
star's  diffuse  environment. 

V.  SUMMARY 

W  e  have  presented  observations  revealing  a  cluster  of  seven 
compact  infrared  sources  in  W40.  all  but  one  of  which  can  Pe 
assiKiated  with  optically  identified  stars.  Being  hcawlv 
obscured  «r.d  centered  on  the  H  u  region,  the  three  hrigntost 
sources  iIRS  la.  2a.  3at  are  reasonable  candidates  for  produc¬ 
ing  ihe  ievel  of  excitation  observed  for  boih  ihe  dust  and  the 
H  ;i  piasma.  The  circumsteilar  regions  observed  for  the  hrignt- 
est  live  sources  i IRS  la.  Ih.  Ic.  2a.  3ai  have  little  impact  on  the 
overail  energetics  id  the  WJi)  dust.  Instead,  the  most  important 
mode  of  excitation  appears  to  involve  the  diffuse  clouds  of 
dust.  The  infrared  energy  distributions  exhiPit  a  wide  ranee  of 
shapes,  suggesting  intrinsic  differences  among  the  distributions 
of  circumsteilar  material. 
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RADIO  AND  INFRARED  OBSERVATIONS  OF  Oil  IR  STARS  AT  THE  TANGENTIAL  POINT 

AND  NEAR  THE  GALACTIC  CENTER 

B.  Bald.1  2  Anm.ila  I.  Sargent,5  M  W.  Wirnir.4  and  A.  F.  Bentley5 

Rtcen  td  IW4  O^tithtr  10.  Jneptej  /V*4  ViHe'm/vr  JO 

ABSTRACT 

We  have  determined  accurate  radio  positions  and  subsequent!)  made  simultaneous  infrared  and  Oil  radio 
measurements  of  a  sample  of  OH  IR  stars  near  the  tangential  point  and  in  the  galactic  center.  Intrinsic  phssi- 
cal  parameters  of  the  stars  and  their  dust  shells,  as  well  as  the  maser  luminosities,  have  been  determined.  Tiiev 
are  consistent  with  the  stars  being  at  the  top  of  the  asymptotic  giant  branch:  the  dense  circumstellar  shells 
suggest  that  they  arc  evolving  rapidly  through  a  phase  of  high  mass  loss. 

Our  sample  comprises  stars  of  widely  different  masses  and  luminosities:  it  is  shown  that  the  OH  luminosity 
is  strongly  dependent  on  the  degree  of  reddening  of  the  circumstellar  shell,  t.e..  on  the  stellar  mass-loss  rate.  In 
addition,  a  quantitative  relation  between  the  OH  luminosity  and  the  optical  depth  in  the  10  /im  silicate  feature 
has  been  determined  for  the  sample.  This  predicts  the  observed  time  variations  in  the  strength  of  the  silicate 
feature  with  time-variable  stellar  luminosity.  While  the  stars  in  the  galactic  center  sample  are  systematically 
less  luminous  than  those  at  the  tangential  point,  the  derived  relatfons  apply  equally  to  each  group. 


Subject  heaJinas:  infrared:  sources  —  masers  —  stars: 

I.  INTRODUCTION 

Combined  infrared  and  OH  observations  are  essential  to 
understanding  the  nature  of  OH  IR  stars,  the  properties  of 
,c  ..  fjj-^jntstellar  envelopes,  and  the  OH  maser  pump  mccha- 
.  tcf.  Werner  er  al  1980:  Herman  19831.  These  cool,  variable 
stars,  typical  of  the  Asymptotic  Giant  Branch  (AGBi  in  the 
H-R  diagram,  have  generally  been  found  through  OH  radio 
surveys  teg,  Johansson  et  at.  1977,  Bowers  1978;  Baud  el  at. 
1979u.  h|,  and  many  are  strong  infrared  sources  with  thick 
circumstellar  dust  shells.  Since  their  OH  maser  emission  can  be 
detected  at  large  distances  from  the  Sun.  OH  IR  stars  can  be 
used  to  study  the  population  distribution,  kinematics,  and 
evolution  of  stars  on  a  galactic  scale  (cf.  Baud  et  at.  1981, 
hereinafter  BHMWt.  Moreover,  these  objects  constitute  a 
unique  tool  for  examining  stellar  evolution  at  the  tip  of  the 
AGB.  a  region  of  the  H-R  diagram  which  is  poorly  understood 
theoretically  and  where  circumstellar  obscuration  often  pre¬ 
cludes  optical  observations.  From  infrared  observations  of 
these  objects,  the  luminosities  of  the  underlying  stars  and  the 
masses  of  the  circumstellar  envelopes  may  be  derived,  w  hile  the 
double-peaked  OH  emission-line  profiles  yield  both  the  stellar 
radial  velocity  r  and  the  circumstellar  shell  expansion  velocity 
c,;  the  latter  quantity  is  statistically  related  to  the  main- 
sequence  mass  iBHMW).  Comparison  of  the  visible,  low  OH 
luminosity  Mira  variables  with  the  optically  obscured,  but 
strong  OH  emitting.  OH  IR  stars  indicates  that  OH  lumin¬ 
osity  is  related  to  the  stellar  mass-loss  rate,  and  it  has  been 
suggested  (Baud  and  Habing  1983.  hereinafter  BH;  Jones  et  at. 
19831  that  such  OH  IR  stars  represent  the  end  point  of  the 
Mira  evolution  on  the  AGB;  their  large  mass-loss  rates  imply 
that  some  mav  be  evolving  toward  a  planetary  nebula  stage. 

In  this  paper  we  present  unambiguous  infrared  identities* 

- .  s)cpt  of  Space  Research.  Gromneen.  The  Netherlands. 

‘  RaJio  Asiron.mu  Lab .  L’lmenitt  of  California.  Berkeley. 

*  Oaens  Valle)  Radio  Observatory.  California  institute  of  Technology. 
Pasadena. 

*  NASA  Ames  Research  Center.  MoiTett  F teld.  California. 

*  Dept  of  Natural  Sciences.  Eastern  Montana  College.  Billings. 


circumstellar  shells  -  stars:  late-type 

tions  of  a  sample  of  OH  IR  stars  in  the  galactic  disk  and  in  the 
_  vicinity  of  the  galactic  center.  Simultaneous  OH  and  broad¬ 
band  infrared  photometric  measurements  allow  us  to  uerive 
quantitative  relations  between  the  infrared  and  OH  duxes  and 
the  properties  of  the  circumstellar  shells,  providing  observa¬ 
tional  evidence  that  these  stars  represent  an  evolutionary 
sequence  of  increasing  mass-loss  rate.  Using  these  relations,  the 
observed  time  variations  of  the  silicate  absorption  feature  with 
changing  bolometnc  luminosity  in  OH  IR  stars  can  be 
explained.  In  addition,  it  appears  that  the  ratio  of  the  OH  flux 
to  the  number  of  pump  photons  is  not  constant  teg..  Ei.uur. 
Goldrcich.  and  Scovilie  1976)  but  is  also  a  function  of  the 
mass-loss  rate. 

ti.  observations 
a)  The  Sample 

A  recurrent  problem  in  determining  the  physical  properties 
of  OH  IR  objects  is  the  uncertainty  and  ambiguity  in  their 
distances,  usually  derived  kinematically  from  the  observed 
radial  velocities.  Most  OH  IR  stars  have  been  discovered  in 
flux-limited  radio  surveys,  so  that  the  near-kinematic  distance 
is  a  plausible  choice  tcf.  Jones.  Hyland,  and  Galley  1983).  but. 
since  the  velocity  dispersions  are  typically  10-30  km  s'1 
(BHMW’l.  the  distance  uncertainties  are  large.  For  example, 
the  near-kinematic  distance  of  OH  36.5  +0.6  is  1.8  kpc.  A 
more  reliable  value  of  I  kpc,  resulting  in  a  sieniiicantly  lower 
estimate  for  the  luminosity,  has  been  derived  using  the 
observed  angular  sire  and  the  light  travel  time  through  the 
shell  (Herman  1983). 

The  present  study  includes  stars  in  the  galactic  plane  whose 
high  radial  velocities  indicate  that  they  are  at  the  tangential 
point— the  position  along  the  line  of  sight  closest  to  the  galac¬ 
tic  center — and  stars  near  the  galactic  center  itself.  For  both 
groups  the  distance  uncertainties  are  relatively  small  The 
tangential  sample  was  selected  from  the  homogeneous  list  of 
BIIMW,  which  is  based  on  systematic  OH  surveys  of  the 
northern  hemisphere  (Johannsson  el  al.  1977;  Bowers  1978: 
Baud  et  at.  !979<i).  It  consists  of  (I)  all  stars  with  radial  velo- 
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-cities  r  larger  than  the  maximum  velocity  allowed  by  cir¬ 
cular  rotation  (cf.  Burton  and  Gordon  197X).  and  (2)  stars  with 
c  >  icmtt  —  10)  km  s  " 1  in  cases  where  v,  >  15  km  s' Stars  in 
the  first  group  arc  assumed  to  be  at  the  tangential  point:  excess 
velocities  presumably  reflect  random  motions.  Although  some 
distance  ambiguity  exists  for  the  second  group,  the  velocity 
dispersion  of  these  stars  is  small  I  -  10  km  s  1 ;  BHMWl,  and 
they  are  probably  within  25  ’.>  of  the  tangential  point. 

The  galactic  center  sample  was  chosen  ( 1 1  from  stars  within  a 
decree  (170  pci  of  the  galactic  center,  with  |rj  >  100  km  s  1 
(Habing  et  al.  19X3)  and  Cl  for  /  >  I  ,  from  a  survey  along  the 
galactic  plane  (Olnon  et  at.  19X1 ;  Baud  el  at.  1979/>|.  The  first 
group  arc  members  of  the  galactic  bulge  population  (e  g.. 
Olnon  ei  at.  I9MI).  where  the  density  of  OH  IR  stars  is  strongly 
enhanced.  The  latter  are  probably  foreground  objects  and 
members  of  the  galactic  disk  population  (see  appendix  A). 

The  distribution  of  t  ,  in  both  samples  is  very  similar  to  that 
for  all  OH  IR  stars  in  the  Galaxy  (eg..  BH.MW).  as  is  the 
distribution  in  galactic  latitude  and  longitude.  However,  since 
their  average  distance  is  large,  these  stars  represent  the  high- 
lummosity  tail  of  the  OH  luminosity  distribution:  L,}H  ranges 
from  1  Jy  kpc2  for  the  nearby  OH-emitting  Miras  to  KXX)  Jy 
kpc:  for  the  brightest  and  most  distant  OH  IR  stars  (cf. 
Nguyen-Q-Rieu  et  at.  1979:  BHi.  but  for  the  samples  here  it  is 
typically  larger  than  100  Jy  kpc:.  Nevertheless,  the  number  of 
stars  in  both  samples  increases  steeply  with  decreasing  OH 
luminosity  (sec  Tables  I  and  2).  which  is  consistent  with  the_ 
OH  luminosity  distribution  found  for  all  OH  IR  stars  (cf. 
BH.MW).  Most  conclusions  drawn  from  the  present  sample 
may  therefore  be  extended  to  the  majority  of  OH  IR  stars. 

b)  Accurate  Radio  Positions  and  OH  Fluxes 

Since  the  OH  objects  in  our  sample  lie  within  a  degree  of  the 
galactic  plane,  and  since  single-dish  OH  measurements  lead  to 


positional  uncertainties  of  z  15'  to  ±  5'.  confusion  with  K  and 
M  giants  along  the  line  of  sight  can  be  a  serious  problem  in 
identifying  their  near-infrared  counterparts  (cf.  Jones  et  at. 
I‘)XI.  19X2i.  Interferometric  observations  of  the  1612  MHz  OH 
maser  emission  were  therefore  carried  out  to  determine  more 
accurate  (±l  "l  radio  positions.  The  need  for  such  positional 
accuracv  is  emphasized  by  the  fact  that  two  stars  common  to 
this  sample  and  to  that  of  Jones.  Hyland,  and  Galley  |19X3) 
appear  to  have  been  misidentihed  by  these  authors  (see  below  i. 

1612  MHz  OH  observations  of  the  tangential  sample  were 
made  on  two  consecutive  days  in  19x0  July  with  the  Very 
Large  Array  (YLA)  spectral  hne  system  of  the  National  Radio 
Astronomy  Observatory.'’  Observations  were  made  with  nine 
telescopes  in  a  hybrid  configuration  and  with  baselines  ranging 
from  0.3  to  3.4  km:  the  synthesized  beam  was  typically 
12'  x  7’.'5.  The  instantaneous  bandwidth  of  780  kHz  was 
divided  among  128  autocorrelation  channels,  each  6.1  kHz 
wide  1 1.1  km  s  '  *).  The  instrumental  phase  and  bandpass  were 
calibrated  using  the  unresolved  continuum  source  1741  -038 
(1  43  Jy  based  on  13.83  Jy  for  3C  286  at  1612  MHzi  and  3C  286 
respectively. 

Positions  were  determined  from  high-resolution  maps  of  the 
high-  and  low-velocity  components  of  the  line  profiles.  In  all 
cases,  emission  from  the  two  components  appeared  unresolved 
and  coincident  within  the  uncertainties  of  the  measurements. 
The  results  are  listed  in  Table  1  in  order  of  increasing  galactic 
loneitude:  column  ll)  contains  the  source  name  and  columns 
(2).  (3).  1 4).  and  |5)  the  positions  and  3  <x  uncertainties  in  R.A. 
and  Deel  respectively.  Although  stars  from  both  this  and  the 
galjcttc-center  sample  are  in  some  cases  displaced  by  0. 1  or 
more  from  the  single-dish  positions,  the  old  source  names  have 

*  NR  AO  iv  operated  bs  Associated  Universities.  Inc.,  under  contract  to  the 
National  Science  Foundation. 
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OH  IR  Suns  AT  THE  Tangential  Point  Am  rate  Radio  Positions.  Velocities,  and  Peak  OH  Fax  Densities 


LV 

HV 

Scs* 

Shv‘ 

1  a 

Name 

R.A  (19501 

3  a 

Deel  (19501 

3  <x 

Ikm  s  ‘ ') 

(km  s"‘l 

(J)l 

(Jyi 

(hi 

in 

(21 

<3i 

(41 

(5| 

161 

(7) 

18) 

(91 

IIOl 

OH  109-1  5  . 

t8"i»r4:'55 

±  ri 

-  18  4(  IS'4 

+  1  2 

i:ou 

144  5 

2.61 

105 

0  (3 

OH  17.2-1.1  . 

18  23  19  55 

+  0  5 

-(4  3008  0 

±0.5 

158  4 

183  3 

1  85 

3  53 

0  13 

OH  182+05  . 

18  19  0721 

+  03 

-12  56  50  2 

±0  9 

113  9 

1377 

1  42 

1  96 

0  12 

OH  13.5  +  1.4  . 

18  16  47.37 

±0.3 

-12  09  27  8 

±0.* 

165  7 

1870 

8  45 

5.50 

020 

OH  20.7  +  01  . 

18  25  44  32 

±0.3 

-10  52  51.1 

+  0  9 

1179 

154  6 

10  8 

6  1 

06 

OH  21  5  +05  . 

18  25  45  49 

+  06 

-10  00  12.4 

+  0.8 

978 

134  0 

25  l 

34  1 

06 

OH  24 .7 -01  . 

18  34  03  6) 

±0  5 

-07  20  52.3 

±11 

109  8 

142  8 

48 

<  1.5 

0  5 

OH  2?  1-0  3  . 

18  35  33.36 

t0  6 

-07  12  34  9 

±14 

1300 

154  8 

7.44 

660 

0  25 

OH  270-0.4  . 

18  39  21  93 

±06 

-05  24  02.8 

±1.1 

87.3 

1164 

3.71 

8  02 

024 

OH  272  +  02  . 

18  37  36.72 

±0.8 

-05  05  28  4 

±0  8 

72.6 

111.4 

2  1 

03 

OH  275-0  9  . 

18  42  01.90 

±06 

-05  12  23.5 

+  0  6 

93  0 

1200 

703 

6  11 

0  25 

OH  27.8  -  (.5  . 

18  44  5797 

±08 

-05  1427.3 

+  0  9 

70  1 

100  5 

2  27 

2.^9 

0  16 

OH  23  5  -  00  . 

18  40  47  45 

±06 

-03  58  57.6 

•  ±0  8 

93  9 

120  6 

11  85 

11  41 

0  24 

OH  29 4-0 8  . 

18  45  12.24 

+  0  4 

-03  32  53  2 

±1.4 

113  6 

1373 

(0  35 

581 

02' 

OH  .10  1-0.7  . 

18  46  04  91 

±0  3 

-02  53  54  t 

±0  8 

78  4 

1190 

57,2 

59  6 

04 

OH  310-0  2  . 

18  46  07.(3 

±08 

-01  51  56.5 

±1  1 

1114 

140  1 

5  22 

6.56 

037 

OH  22.0-0  5  . 

18  48  51.21 

±0.8 

-01  07  29.3 

+  0  6 

550 

96  0 

83 

51 

04 

OH  J2  1+09  . 

18  44  04  62 

+  08 

-00  20  29  9 

+  0  6 

124  7 

149  3 

2.36 

3.77 

017 

0(1  34  9  +  08  . 

18  49  43  87 

±1.1 

+  02  0008  0 

+  0  7 

540 

82.9 

2  39 

3  (0 

Oil 

OH  36  4  +0  3  A  . 

18  54  Ot.52 

+  0  6 

+  03  1608.2 

±0  6 

83  5 

121. 1 

4  88 

2.76 

0  19 

OH  36  4+0  3  B  . 

18  53  5761 

+  03  1321.1 

OH  364 +03 C . 

18  54  0715 

+03  1402.1 

OH  371  -08  . 

18  59  36  22 

±08 

+  03  15  53.3 

±11 

75  0 

101  5 

1441 

20  32 

047 

OH  37.7-1.4  . 

19  02  4007 

±0  3 

+  03  36  23.4 

+  0.6 

tooo 

1199 

2.32 

0.47 

0  16 

OH  39.9-0.0  . 

19  01  42.93 

±06 

+  06  08  45  0 

+  0  8 

1337 

163 1) 

3  6 

6.2 

04 

•  1612  OH  Dux  density  at  the  limeo(  Ihe  infrared  observations:  spectral  resolution  is  A  9  kHz. 
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)  TABLL  2 

■A.  ’ 

OH  IR  Stars  in  niE  Galactic  Ctntir  Region: 

Arc  innR  cum  Pum  thins.  Veldt  mis.  and  Peak  OH  I  u  s  Oi  nsitiis 


Name 

III 

R  A  119501* 
(2l 

Dccl  ||950|* 
131 

LV 

(km  v  ') 
(4) 

HV 

(km  s  ') 
(5l 

S,v* 

(J>l 

161 

s  * 

Uvl 

P) 

OH  359  1  +  1  1 . 

I7*35'56-9X 

-  29  02  24  8 

-146  5 

-128 

2.3 

2.8 

OH  359  2  +  0.2 . 

17  39  55.33 

-29  29  36  1 

-150 

-121 

09 

20 

OH  359  4  +  0  1  . 

17  40  34  13 

-29  24  59  5 

- \w 

-223  5 

08 

22 

OH  02 +  00  . 

17  42  45  49 

-28  44  10  0 

+  145 

+  174 

1  t 

07 

OH  0.3 -0.2  . 

17  43  56  64 

-28  43  41  2 

-356 

-327 

2.3 

19 

OH  0  5 -0.2  . 

17  44  14  44 

-28  35  3!  8 

+  130 

+  153  5 

09 

1.5 

OH  1  OS +  0.4 . 

17  43  35.37 

-27  48  47  2 

-142 

-  108 

<06 

<06 

OH  1.1-0*  . 

17  48  16 .88 

-28  24  52  7 

-9 

+  30 

20.6 

24  5 

OH  15-0 1  . 

17  46  12  10 

-27  41  O0  8 

-141  5 

-114  5 

27 

5  1 

OH  1.7 -0.0  . 

17  46  53.55 

-27  24  10  0 

-105 

-134 

0.8 

3.5 

OH  26-04  . 

17  50  11.14 

-26  5601  5 

-27 

+  18  5 

13  3 

14  2 

OH  40  +  09  . 

17  48  17  90 

-25  01  08.7 

+  59 

+  88 

1.3 

10 

OH  50-04  . 

17  54  32  21 

-25  12  42.6 

-  100 

-  146 

1.8 

25 

OH  90-0.1  . 

18  0  1  09  13 

-21  14  02  4 

—  64 

-33 

<0.9 

<09 

OH  9.6  -0.5  . 

18  02  10.19 

-  20  22  30  6 

—  49 

-76 

20 

35 

OH  12.8-09 . 

18  13  53  45 

-18  16  08  9 

-67 

-44 

97 

2.3 

OH  12.8  +  0.3 . 

18  09  23  U9 

-  17  4348  I 

+  125 

+  152 

<  10 

<  10 

OH  14  8+0.2 . 

18  13  48  70 

-16  05  15  0 

+.109 

+ 133 

<  1  0 

<10 

*  3  a  uncertainties  in  R  A.  »  jO'3  and  in  Ded.  =  s0'9 

*  1612  MHz  OH  rtus  density  on  1981  Aug  2?-Sep  I:  upper  limits  ire  3  o.  spectral  resolution  is  49 
kHz.  3  c  unccrtJiniv  is  r^2J>. 


been  adopted  for  consistency  with  other  published  work. 
Limited  LV  coverage  for  OH  36.4  *  0.3  resulted  in  some  uncer¬ 
tainty:  position  A  is  the  most  likely  but  two  less  plausible 
.  ssitions  (B  and  Cl  are  also  given.  The  positions  of  OH 
' i'.l—  1.4.  observed  on  both  days,  agreed  within  072,  confirm¬ 
ing  the  accuracy  of  the  measurements.  Comparison  w  ith  more 
recent  VLA  observations  of  some  of  the  stars  (Herman  1 983k 
further  strengthens  our  confidence  in  the  results.  For  OH  20.7 
and  OH  31.0  the  near-infrared  positions  listed  by  Jones. 
Hyland,  and  Gatley  (19831  arc  more  than  10"  offset  from  our 
improved  radio  positions,  suggesting  that  their  infrared  objects 
may  not  be  associated  with  the  OH  masers. 

Accurate  positions  for  the  galactic-center  sample  were  deter¬ 
mined  with  the  VLA  on  1981  May  7.  At  this  time,  by  using  a 
bandpass  of  390  kHz  and  64  autocorrelation  channels  of  6  1 
kHz  each.  13  teiescopes  could  be  employed,  covering  a  range  in 
baselines  from  0.2  to  10.5  km.  The  synthesized  beam  was 
T  x  4".  Instrumental  phase  calibration  was  effected  using  the 
standard  calibration  source  1748-253  and  the  low -elevation, 
secondary  calibration  source  1741  +038.  The  names  of  the 
stars  and  their  radio  positions  are  listed  in  Table  2  in  columns 
(1).  (2).  and  (3).  The  data  on  1741  —038  show  that  the  positional 
uncertainty  due  to  systematic  phase  errors  is  less  than  071. 

Single-dish  OH  (luxes  of  the  tangential  stars  were  measured 
with  the  Dwingcloo  25  m  radio  telescope  in  1980.  within  20-60 
davs  of  the  Wyoming  infrared  observations  (see  tj  lie).  Since 
this  time  difference  is  much  less  than  10%  of  the  ty  pical  periods 
for  these  stars  (Herman  1983),  the  infrared  and  radio  fluxes  can 
be  considered  coeval.  Fluxes  for  the  galactic-center  sample 
were  measured  with  the  85  foot  (26  ml  telescope  at  Hat  Creek 
between  1981  August  27  and  September  I.  one  month  after  the 
infrared  observations  on  the  IRTF  (see  ij  IK).  No  single-dish 
I  observations  were  made  at  the  lime  of  the  UKIRT  infrared 
ooserving  run  in  1982. 

The  high  spectral  resolution  of  the  single-dish  measurements 
(0  46  km  s  ~ ')  provided  very  accurate  velocities  and  peak  lluxcs. 
Tables  I  and  2  contain  the  measured  radio  parameters  of  the 


-  two  velocity  components  of  the  OH  emission  profiles  for  the 
tangential  and  the  galactic-center  samples  respectively:  radial 
velocities  with  respect  to  the  local  standard  of  rest,  of  the  low- 
and  high-velocity  (LV.  HV)  peaks  of  the  OH  line  protile  are 
presented  in  columns  (6)  and  (7i  of  Table  1  and  columns  |4|  and 
(5i  of  Table  2;  the  corresponding  peak  flux  densities.  i‘LV  and 
SHV.are  shown  in  columns  iS)  and  (9)  of  Table  1  and  columns 
(6i  and  (7|  of  Table  2.  The  1  a  uncertainty  in  the  OH  flux 
density  measurement  of  the  tangential  sample  is  listed  in  the 
last  column  of  Table  1. 

c)  Infrared  Observations 

Infrared  observations  of  the  tangential  sources  were  made  at 
the  Cassegrain  focus  of  the  Wyoming  Infrared  Telescope 
(WIROi  in  August.  19"'9  and  1980.  using  a  heiium-cooled 
bolometer  iBentiey  1980).  The  aperture  and  chopper  throw 
were  775  and  20”  respectively.  Infrared  observations  of  the 
galactic-center  sample  were  made  at  the  IRTF  and  UKIRT  on 
Mauna  Kea.  Hawaii,  in  1981  July  and  1982  June  respectively, 
using  the  standard  helium-cooled  bolometers  and  also  the 
liquid  nitrogen-cooled  InSb  photometer  at  the  UKIRT.  A  6" 
aperture  was  employed  at  the  IRTF  and  one  of  775  at  the 
UKIRT:  chopper  throws  were  15"  EW  at  IRTF  and  30"  in 
declination  at  UKIRT.  Time  did  not  allow  us  to  observe  all  the 
stars  in  Table  2. 

In  general,  measurements  of  the  infrared  counterparts  of  the 
OH  IR  stars  were  made  by  integrating  on  the  VLA  positions 
using  a  broad-band  filter  centered  at  10  pm.  where  the  energy 
distributions  usually  peak  and  where  confusion  by  foreground 
K  and  M  giants  is  negligible.  However,  for  a  few  galactic -center 
sources,  searches  at  UKIRT  were  carried  out  at  8.7  or  20  pm 
when  the  earlier  IRTF  observations  had  shown  that  either  the 
10  pm  absorption  feature  was  sufficiently  strong  to  preclude  a 
10  pm  detection  or  that  the  energy  distribution  peaked  long- 
ward  of  12.5  pm.  Infrared  positions  determined  for  the  stronger 
sources  from  a  five-point  grid  centered  on  the  VLA  positions 
agreed  well  with  the  radio  positions.  Long  integrations  adjacent 
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W|RO 

Zeroth  Mao 
f, 

(Jvl 

IR  TF 

UkIRT 

Ziroth  Mao* 
f. 

(3)1 

A 

A/. 

Ipmi 

A 

|/jrm 

A/. 

(pmi 

A 

1*4  ITI! 

A/. 

Ipml 

23 

0.7 

59.S 

36 

t  3 

277 

38 

067 

3.8 

065 

238 

49 

0  7 

15S 

48 

0.57 

4.X 

060 

153 

87 

10 

54 

8.7 

1.19 

8.7 

12 

50 

9  R 

10 

9.7 

1.0 

40 

10  0 

5.x 

41 

10  3 

103 

10  3 

1.25 

34 

114 

20 

32 

116 

1.26 

116 

1.17 

28 

136 

08 

26 

125 

1.17 

125 

125 

24 

19  5 

58 

12 

20 

90 

20.9 

60 

10 

*  Appii«  lo  both  LMRT  and  IRTF. 


to  the  radio  positions  of  a  feu-  objects  yielded  no  detections. 
We  are  therefore  confident  that  the  infrared  sources  can  be 
identified  unambiguously  with  the  OH  stars. 

Following  the  detections  of  the  IR  counterparts,  photo¬ 
metric  measurements  between  3  and  20  were  made.  Appro¬ 
priate  filter  characteristics  are  presented  in  Table  3.  After 
taking  account  of  the  relevant  air-mass  corrections,  magni¬ 
tudes  at  each  wavelength  were  determined  from  comparison  - 
with  standard  stars.  Conversion  from  magnitudes  to  flux  den¬ 
sities  was  accomplished  using  the  values  shown  in  Table  3  icf. 
lehrz.  Hackwell.  and  Jones  1974.  Beckwith  ef  a/.  1976i.  The 
measured  infrared  flux  densities  and  3  a  upper  limits  for  the 
tangential  sample  are  given  in  Tables  4  and  5.  while  those  for 
the  galactic-center  sample  are  in  Table  6.  In  both  tables,  the 
year  in  which  the  observations  were  made  is  indicated  in 
column  (2).  Corresponding  infrared  energy  distributions  are 
sh'swn  in  Figures  I  and  2  respectively.  For  wavelengths 
between  3.8  and  12.3  rim  the  uncertainties  in  the  fluxes  are  of 
order  z  10%.  and  at  20  itm  about  15%.  Typical  error  bars  for 
the  shorter  wavelengths  are  plotted  at  the  8.7  jim  data  point  in 
each  source  energy  distribution.  Error  bars  for  the  20  #im  flux 
are  also  indicated. 


Upper  Lours  fur  Tvm.intial  Stars 
(1*1X0  Observations! 


10  (.m  Flux  Density 


Name  (Jyl 


OH  1?: .  <39 

Oil  18.5 .  <3.5 

OH  35  1 .  <35 

OH  3%) .  <33 

011  37  5  .  <4  5 

01(39  4 .  <40 

OH  lb 4 .  <3.6 

OH  37  | .  <3.1 


tit  RESULTS  AND  ANALYSIS 

Long-term  monitoring  observations  at  18  cm  (Harvey  et  al. 
1974;  Jewell.  Webber,  and  Snyder  1980:  Herman  1983)  show 
that  L,)h  and  L,  change  simultaneously,  strongly  suggesting 
that  the  OH  maser  is  radiatively  pumped.  But  from  the 
observed  correlation  between  LOH  and  the  shell  size,  and  from 
the  very  broad  OH  luminosity  distribution  (cf.  BH.MWl,  BH 
have  argued  that  L,)H  is  a  steep  function  of  the  stellar  mass-loss 
rate  and  more  or  less  independent  of  L.  This  is  also  suggested 
by  the  observed  correlation  between  LOH  and  near-infrared 
colors  (eg..  Jones.  Hyland,  and  Gatley  1983;  Engeis  1982: 
Herman  1983).  This  raises  the  question  of  the  relation  between 
M  and  Loh.  and  the  nature  of  the  underly  ing  star.  As  a  result  of 
(1)  homogeneity.  (2)  simultaneous  radio  and  infrared  observa¬ 
tions.  and  (3)  information  at  /.  >  10  ;im.  the  present  sample 
allows  a  quantitative  analysis  of  the  data  involving  LOH.  L,, 
and  M  and  a  better  understanding  of  the  influence  of  M  on  the 
maser  pump  efficiency.  The  physical  parameters  of  the  stars 
and  their  circumsteilar  shells,  determined  from  the  radio  and 
IR  observations  described  above,  are  summarized  in  Tables  7 
and  8.  Their  derivation  and  the  influence  of  extinction  are 
discussed  in  Appendices  A  and  B  respectively. 

o)  The  OH  Maser  Flux  and  Stellar  Luminosity 
L0h  ranges  from  1-10  Jy  kpc:  for  nearby,  optically  identified 
OH  emitting  Miras  to  10-1000  Jy  kpcJ  for  the  distant.  OH 


TABLE  4 

Infrared  Flux  Densities  of  Tangential  OH  IR  Stars 
(Jvl 


Wavelength  i«mi 


Name 

Date 

2.3 

36 

49 

8.7 

10.0 

114 

12  6 

195 

OH  10  8  ♦  1.5 . 

1980 

01 

11.7 

38 

39 

129 

29.1 

OH  182  +  0.5 . 

19X0 

0  7 

107 

8.1 

6  1 

152 

<12  8 

OH  20.7  *•  0  1 . 

1979 

<0.2* 

08 

2.3 

2.5 

0.9 

5.8 

7.9 

OH  21  5+0  5 . 

1979 

<01* 

2.1 

206 

350 

29  2 

16  7 

634 

68.5 

OH  24  7-0  1  . 

1979 

09 

116 

27.7 

305 

24  5 

22.0 

336 

28  3 

OH  27  2+0  3 . 

19-9 

02 

2.2 

50 

83 

78 

8.8 

10.7 

92 

OH  278-1.5 . 

1980 

32 

147 

38  8 

37.7 

37.2 

38  6 

316 

OH  3X5-00 . 

19X0 

03 

6.5 

7.4 

<6  8* 

OH  301-0  7 . 

1979 

<01* 

14 

169 

68  6 

53.6 

28.7 

100.5 

129  4 

OH  <10-0.2 . 

19X0 

<0.2* 

19.8 

14  1 

<14  7* 

24  8 

194 

OH  32.0-0.5 . 

1979 

<01* 

0.1 

2.8 

5.9 

64 

39 

136 

15.6 

OH  32.1+09 . 

1980 

01 

1.5 

120 

<7.0* 

11.4 

OH  34  9  +  0  8 . 

19X0 

<01* 

1.8 

101 

66 

10.3 

OH  37  7-  1  4 . 

19X0 

<01* 

II 

37 

7.1 

83 

<5  4* 

OH  399-00 . 

1979 

<0  1* 

l.l 

49 

116 

92 

68 

154 

15  6 

•  Upper  limits  are  3  a 


63: 

) 


W  iHUV.tii  i/.mi 


Nvwf 

Dvir 

J.X 

4.9 

8  7 

97 

105 

115 

125 

19  H 

011  159 .1  ♦  1.1 . 

i9x: 

A 

C 

ia 

008 

0.31 

014 

104 

3.13 

19X1 

Oil 

0  61 

0.73 

0  79 

161 

3  31 

10  67 

011  <59  4  *  0  | . 

|9X2 

<o:* 

0  45 

<16* 

19M 

o:3 

o.;: 

0  37 

0  56 

1  87 

331 

oho: +00  . 

i9x: 

095 

1  24 

1  73 

0.18 

031 

1.73 

1  85 

2  60 

|9S| 

0.58 

047 

0.0*1 

<0  2* 

060 

0.5 

0  50 

OH  0  3  -0.2  . 

WSI 

0  45 

12 

OH  0  5-02  . 

1 9*2 

033 

070 

<04* 

073 

1.10 

1.27 

19X1 

051 

104 

o:4 

029 

1.18 

3.11 

32 

011  1  OS +0.4 . 

19X1 

<0  3* 

<0  3* 

023 

42 

on: p -04  . 

i9x: 

78  09 

80  30 

197  33 

159  34 

167  39 

19194 

113.78 

167  49 

OH  40-0.9  . 

19X2 

<0.2* 

0  63 

0.84 

<2.1* 

5.25 

OH  4.6  —0.4  . 

i9x: 

360 

493 

853 

1.36 

1  63 

7.78 

9  91 

8  95 

I9>l 

644 

806 

2  11 

233 

923 

17  42 

23.3 

0119  6  +  0  5  . 

i9x: 

7.88 

644 

i8.i: 

14  91 

15  40 

11.83 

1543 

1585 

19X1 

593 

167’ 

19  33 

17  88 

19.52 

19  OX 

170 

*  L.  pper  limns  are  3  <r. 
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luminous,  but  optically  obscured.  OH  IR  stars  studied  in  this 
paper.  The  OH  luminosity  of  these  objects  increases  steeply 
with  decreasing  i.OH.  and  it  has  been  suggested  that  Miras  may 
-■present  the  low-luminosity  tail  iBHMW'i.  Values  of  L,  for 
liras  and  OH  IR  stars  are  not  significantly  different,  leading 
to  the  hypothesis  that  L,,H  is  independent  of  L.  (BH)  Indeed, 
the  present  data  support  this  contention:  in  Figure  3.  where  the 
tangential  sources  are  plotted  as  tilled  circles  and  the  galactic- 
center  sampie  is  represented  by  crosses,  there  is  no  evidence  of 
a  strong  correlation  between  LOH  and  Lm. 

E„„  for  the  galactic-center  sample  is  a  factor  of  2  to  3  smaller 
than  that  for  the  tangential  sample,  reflecting  the  higher  sensi¬ 
tivity  of  the  OH  observations  near  the  galactic  center.  The 
mean  bolometric  luminosity  for  the  detected  stars  in  the 
tangential  sample  is  3.4  x  I04  L  -,.  whereas  it  is  2  x  101  Ls  for 
the  stars  in  the  vicinity  of  the  galactic  center  (/  <  l  I.  It  is  not 
clear  whether  this  discrepancy  is  real,  ft  could  result  from  dif¬ 
ferences  in  the  sensitivity  and  completeness  of  each  sample  or 
from  the  variations  between  the  space  density  of  OH  IR  stars 
near  the  galactic  center  and  in  the  disk  near  tangential  point. 
Although  OH  IR  stars  in  the  galactic  center  have  much  lower 
bolomeiric  luminosities  than  those  at  the  tangential  point,  a 
larger  interstellar  extinction  toward  the  galactic  center  could  to 
some  extent  explain  the  disparity  .  Ft  is.  however,  unlikely  that  it 
can  account  for  a  discrepancy  of  more  than  a  factor  of  2. 
indicating  that  the  OH  IR  stars  in  the  bulge  are  intrinsically  at 
least  5  times  less  luminous  than  those  in  the  disk. 

Figure  4  demonstrates  that  there  is  no  correlation  between 
L„„  and  t  ,.  the  shell  expansion  velocity  .  Given  the  absence  of  a 
correlation  in  Figure  3.  this  is  scarcely  surprising  because,  for 
OH  IR  stars,  i,  is  statistically  related  to  the  main-sequence 
mass  lUHMW  |.  Since  the  main-sequence  mass  is  related  to  the 
4GB  mass  (of.  Ibcn  1981k  which  in  turn  determines  the  bolo- 
.■tric  luminosity  (Paczynski  1970).  r,  and  Lm  should  in  fact  be 
correlated.  The  statistical  data  of  Jones.  Hyland,  and  Galley 
1 198  J)  appear  to  support  such  a  correlation  and.  although  the 
scatter  in  Figure  5  is  large,  there  is  a  weak  trend  for  the  most 
luminous  stars  to  have  higher  values  off,. 


b  I  T/ieOH  M  user  Flux  and  I  he  Circumstellur  Dust  Density 

The  long-term  variability  and  bolomeiric  luminosity  of  these 
OH  IR  stars  suggest  that  their  stellar  parameters  are  sinvlar  to 
those  of  the  OH  emitting  Miras  with  effective  temperatures 
between  2500  and  3500  k  and  maxima  in  the  energy  distribu¬ 
tions  between  I  and  2  pm  tef  Merriil  1977|.  It  is  commonly 
thought  that  the  OH  maser  is  pumped  by  35  pm  photons 
(Elitzur.  Goldreich.  and  Scoville  1976).  but  in  the  case  of  the 
OH  luminous  OH  IR  stars,  the  stellar  photospheric  flux  at  35 
pm  is  clearly  insufficient  to  provide  enough  maser  pump 
photons.  This  implies  that  the  maser  luminosity  depends  on 
the  efficiency  with  w  hich  the  stellar  light  is  converted  to  35  urn 
emission  in  the  circumstellar  shells,  i.e..  on  the  dust-column 
density  If  the  pumping  is  due  to  35  pm  photons  and  the  maser 
is  saturated,  then  stars  with  the  same  OH  luminosity  but  differ¬ 
ent  values  of  L.  produce  similar  35  pm  fluxes,  and  the  energy- 
distributions  must  be  redder  and  cooler  for  lower  values  of  L%. 
Thus  the  ratio  LOH  L,  should  be  positively  correlated  with  the 
fraction  of  the  total  energy  that  is  emitted  ai  35  pm. 

In  Figure  6  Lim  L»  is  plotted  against  the  fraction  of  the  total 
energy  emitted  at  20  pm.  /.F» 0  F,,.  Since  in  most  cases  the 
maximum  in  the  energy  distribution  of  these  objects  lies  short- 
ward  of  35  pm,  >.F,0,  F,«  is  a  reasonable  measure  of  the  relative 
amount  of  energy  emitted  at  35  pm.  Both  LOH  L .  and  /.F,0  F,* 
arc  distance-independent.  The  OH  luminosity  increases  steeply 
with  the  fraction  of  the  total  emission  that  emerges  at  20  pm. 
that  is.  w  hen  the  energy  distribution  becomes  redder.  It  is  clear 
from  this  figure  that  for  most  stars  there  is  a  remarkably  tight 
correlation  between  the  two  quantities,  which  is  intrinsic  and 
cannot  be  due  to  variable  interstellar  extinction.  Three  objects, 
OH  359.1,  OH  4.6.  and  OH  10.9.  which  lie  significantly  below 
the  observed  correlation,  have  anomalously  low  OH  lumin¬ 
osities.  OH  0.2  lies  well  above  the  trend  and  appears  to  emit 
more  OH  photons  than  35  pm  pump  photons,  in  violation  of 
the  simple  35  pm  pump  model  I  Elitzur.  Goldreich.  and  Scoville 
1976).  This  discrepancy  is  discussed  in  IV. 

A  linear  regression  fit,  excluding  the  four  anomalous  sources 
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ha.  2. — Infrared  energy  distributions  of  ihe  galactic  center  sample  measured  *nh  l RTF  m  19X1  and  *ifh  L  KIRT  in  19X2  iJjsheJ) 


with  F,„  and  /F.0  in  units  of  10' 13  W  m  ' 2.  The  dependence  of 
SOH  on  Flt  and  /.F:o  F,t  is  consistent  with  the  observations; 
the  OH  maser  output  follows  closely  the  time  variations  in  L.. 
as  would  be  expected  from  SOH  x  F,*  In  addition,  the  steep 
dependence  of  SOH  on  /.F;oFi*  implies  that  no  direct  correla¬ 
tion  is  expected  between  L,  and  tcf.  Fig.  3l:  the  relative 
amount  of  long-wavelength  flux,  which  is  a  function  of  the 
dust-column  density  and  hence  the  mass-loss  rate,  rather  than 

TABLE  7 


Raisio  and  Ivr*A«FO  Pmipfunrs  or  Tavjntiai  OH  IR  Stars 


Name 

Distance' 

Opel 

Us*. 

(Jy  8 pc'! 

(km  s'1) 

F  m 

( 10 '  ‘  *  V.  m  "  *1 

f; 

(10*  L,.) 

T, 

1 K.  > 

'..i 

OH  10.9  ♦  15 . 

85-118 

120-230 

16  25 

63  £  7 

1.34-258 

310 

1.3 

OH  172-1.1 . 

96 

240 

12.50 

OH  18  2*05 . 

95 

150 

1215 

50  ±  6 

1.33 

410 

09 

OH  185*14  . 

95 

620 

1065 

OH  20  7*01 . 

94 

720 

18  30 

18  +  2 

A  47 

350 

18 

OH  21  5*05 . 

8  6-101 

2pd-:990 

18  20 

260  ;  30 

5  68-7  83 

410 

1.2 

011  24  7-0  1  . 

9.1 

220 

16  50 

265  »  30 

6.48 

580 

04 

OH  25  1  -03 . 

9  1 

590 

1240 

OH  27  0  —  0.4 . 

7  4-10  4 

30l  >-590 

1460 

OH  272*02 . 

6  7-11.1 

90-260 

19  40 

60  £  7 

080-218 

520 

0.2 

OH  275-09 . 

89 

520 

13  50 

OH  278-  1  5 . 

6  1-116 

90-340 

15  20 

340  t  40 

3.73-13.5 

620 

<0  1 

OH  28  5  -  00 . 

88 

900 

13  35 

49  +  6 

1.12 

•>380 

OH  294  -  08  . 

8  7 

590 

11.85 

OH  30  1  -0  7 . 

7.3-IOO 

3110-5840 

20  30 

380  *  45 

5  98-11.2 

360 

1.2 

OH  310-02 . 

86 

430 

14  35 

96  -  10 

2.11) 

<340 

>04 

OH  120-05 . 

54-116 

190-880 

20  50 

45  i  5 

0  39-1  79 

410 

10 

OH  321  *0.9 . 

8  5 

220 

12.30 

75  *  9 

1.60 

500 

OH  34  9  *  08 . 

-*8-11.6 

60-370 

14  45 

56  £  7 

0  38  2.22 

490 

0.5 

OH  36  4  +  0.3 . 

8.1 

240 

18  85 

OH  37  1  -  08 . 

80 

240 

1325 

OH  37  7-1  4 . 

79 

70 

1000 

56  £  7 

1  03 

>490 

OH  399-00 . 

7  7 

280 

1465 

70  £  8 

1  23 

430 

08 

*  In  case  of  a  distance  ambiguity,  values  for  both  the  near-  and  the  far-lunematic  distance  are  given. 


•nd  the  upper  limit  for  OH  18  2.  yields 

log  -  o  n  log  -  o.5t±o.o3> .  <d 

The  harmonic  mean  OH  flux.  ^OM*  is  then  given  by 

/;  f  V* 

5oh  =  0  Jv  •  <-> 


T 
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Radio  and  Inharid  Pnoprants  or  nit  GALAf-nr  Cinttr  OH  IR  Sr»*s 


Name 

Distance 

Ik  pet 

<J>  kpcj) 

(km  s  ‘  *1 

F„ 

HO1*  W  m  ’) 

L. 

(10*  L  I 

Tt 

(Kt 

Near 

OH  .'59  1  ... 

to 

925 

54 

0.16 

340 

|4k  : 

250 

126 

037 

300 

02 

I9X| 

OH  <594  ... 

to 

1225 

<  14 

<004 

<230 

|Yx2 

130 

<4.5 

<0  13 

<370 

>10 

19X1 

oho: . 

to 

14  5 

119 

0  34 

590 

!  t) 

mz 

90 

3.4 

0  1 

450 

2.1 

19X1 

OH  O'  . 

10 

210 

14  5 

>1  2 

>0  *)4 

19X1 

011  0  5  . 

to 

11.75 

38 

0.11 

460 

>08 

19x2 

120 

8  1 

0  24 

400 

>  1.7 

19X1 

OH  I  08 . 

10 

<  60 

170 

<48 

<0  14 

19X1 

OH  2  6  . 

5 

340* 

22  75 

936 

68 

600 

02 

|9.\2 

OH  4  0  . 

to 

HO* 

14  5 

<  6  6 

<0  19 

370 

1982 

OH  4  6  . 

10 

230 

56 

1  6 

530 

1  8 

19S2 

210 

89 

26 

500 

16 

19X1 

OH  9  6  . 

to 

23.5 

90 

26 

50t) 

02 

mz 

260 

98 

29 

461) 

<0  1 

19X1 

*  I9X|  measurement. 


the  integrated  flux  (stellar  luminosity),  is  the  dominant  par¬ 
ameter  governing  the  maser  flux.  That  is  to  say.  the  OH  lumin¬ 
osity  depends  mainly  on  the  effective  temperature  of  the 
circumstellar  dust  (see  $  Uh  l.  with  an  overall  sealing  factor  due 
to  the  bolometric  luminosity.  Since  the  precise  relation 
between  the  relative  amount  of  20  jim  flux  and  M  is  not 
known,  a  direct  comparison  between  LotlL,  and  M  is  not 
possible. 


c)  Silicate  Absorption 

The  silicate  absorption  feature  at  9  7  um  is  formed  in  the 
cool  outer  regions  of  the  circumstellar  shell  Its  apparent  depth. 
t„,.  is  a  measure  of  the  dust-column  density  and  is  also 
expected  to  be  inversely  proportional  to  T.  ikwan  and  Scoville 
19761.  Since  r,,,  appears  to  be  correlated  with  Loli  L.  (Hetman 
19S3I.  it  should  increase  with  /.F:0  F,„:  Figure  la  shows  that 
this  is  indeed  the  case,  although  OH  0  2  and  OH  4.6  have 
unexpectedly  high  values  of  r„,.  Both  objects  are  also  anom¬ 
alous  in  Figure  5.  OH  359  1  is  extremely  red  but  has  a  low 
value  of  t„,.  similar  to  OH  17.7-2.0  lOlnon  el  al.  1 98-41.  sug¬ 
gesting  that  in  the  coolest  objects,  which  presumably  experi¬ 
ence  the  largest  mass-loss  rate,  the  formatior  cf  the  silicate 
feature  is  severely  inhibited  by  an  as  yet  unspecified  mecha- 


Fio  3  —The  stellar  bolometnc  luminosity  L,  in  L~  as  a  function  of  Ihe  maser  lummosuv  L rm  in  Jv  kpcJ  In  ihis  and  the  followm*  figures  ihe  unccnnal  sample  is 
indicated  by  tilled  circles.  Ihe  galactic  cenier  sample  by  crosses  Vjlues  corresponding  to  ihe  near-  and  (jr-lunemjtic  distances  lor  ihe  tangential  -ample  are 
connected  by  a  line  The  arrow  imi-cjies  ihe  sensitivity  limit  for  Ihe  W(RC)  measurements  isee  tes! I 
hio  4  — The  maser  luminosity  as  s  function  of  i,.!he  outflow  yeloctly  of  the  circumstellar  shell 
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:9. 


log  (XF2S/Hr) 

Fiu.  6 


F Hi.  5  —  B.ilometric  luminosity  as  a  function  of  outflow  velocity  — 

F:i.  6  —  The  ratio  ol  nu<r  (o  buiometnc  luminosity  j>  a  function  of  the  traction  of  tout  infrared  energy  emitted  jt  20  rnn  The  open  circle  represents  OH  I'll  .ec 
•evil  L.,h  and  L ,  are  in  the  same  units  as  in  F  ig  .1  Inirured  energies  are  in  Wm  ■ :  The  dasned  a'row  indicates  tnc  effect  of  2?  mag  oi  sisuai  extinction  isee  Appendix 
The  soud  line  is  a  least-sijuares  tit  through  the  tangential  sample,  excluding  the  anomalous  sources  isce  eq.  It  Typical  error  bars  are  indicated 


nism  A  linear  tit  through  all  data  points  except  the  limits  and 
the  two  anomalous  stars  gives  the  following  relation: 

2.11*0.2)  log  +  l.SrO.2).  (3) 

with  a  correlation  coefficient  of  92"«.  Figure  7b  confirms  that 
the  silicate  optical  depth  is  proportional  to 

An  tnsterstellar  component  in  the  silicate  absorption  feature 
has  not  vet  been  completely  ruled  out  In  a  small  sample  of 
OH  I R  stars  Evans  and  Beckwith  ( 1977*  found  no  correlation 


between  r„,  and  distance,  arguing  for  a  circumstellar  origin.  If 
the  silicate  feature  were  due  solely  to  interstellar  extinction,  no 
correlation  with  /.F:o  F„  would  be  seen  in  Figure  7a  and. 
indeed.  would  be  expected  to  depend  on  galactic  latitude. 
Since  this  is  not  the  case,  the  possibility  of  the  silicate  feature 
arising  in  the  interstellar  medium  can  be  ruled  out. 

IV  DISCUSSION 

Our  simultaneous  radio  and  infrared  observations  of  two 
samples  of  distant  OH  IR  stars  at  the  tangential  point  and  the 


l°g(XF20/FIR)  TC(KI 

Ftti.  7  —The  apparent  silicate  optical  depth  as  a  function  of  lulthe  fraction  of  infrared  energy  emitted  at  CO  «m  and  iM  the  color  temperature  between  4  8  and  12.5 
/im  The  solid  line  m  lui  is  a  least-squares  tit  lo  the  tangential  sample  (see  eq  Jl 
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galactic  center  demonstrate  that  the  OH  luminosities  of  these 
objects  are  dependent  on  the  fraction  of  the  total  infrared 
energy  which  is  emitted  near  35  jim.  which  in  turn  is  a  measure 
of  the  mass-loss  rate.  Sources  with  the  same  stellar  parameters, 
such  as  bolometnc  luminosity  and  mass,  may  have  widely  dif¬ 
ferent  circumstellar  properties.  This  is  reflected  in  large  varia¬ 
tions  in  Loi,  among  stars  with  similar  values  of  t,  and  For 
a  given  bolometnc  luminosity,  stars  with  higher  Oil  lumin¬ 
osity  have  redder  infrared  energy  distributions  and  deeper  sili¬ 
cate  absorption  features,  indicative  of  cooler  circumstellar 
shells  and  higher  mass-loss  rates. 


pump  photons  arc  produced  grows,  and  ,ir})  increases  due  to 
the  changing  geometry  between  the  maser-cmitting  region  and 
the  35  jim  pump  region  tcf.  Evans  and  Beckwith  1977). 

hi  I'ahutwns  in  the  Silicate  Absorption  Depth 
Our  results  also  indicate  that,  for  a  given  star,  as  L,  declines, 
an  increase  in  the  silicate  absorption  feature  r,„.  should  result. 
W  hile  the  luminosity  of  the  OH  maser  changes  in  phase  with 
the  bolometnc  luminosity  of  the  central  star,  the  amplitude  of 
the  variation  is  about  a  factor  of  2  smaller  (e  g..  Herman  1983). 
with 


a)  Maser  Pump  Efficiency 

Equations  (I)  and  (2l  show  that  for  a  given  star  L,)H  is  a 
function  of  the  relative  amount  of  long-wavelength  photons 
produced  in  the  circumstellar  shell.  In  the  35  gm  pumping 
scheme  the  observed  maser  pump  efficiency  t  is  given  by 

$OM  , 

*  -  c  /'15  r 

(cf.  Evans  and  Beckwith  1977).  where  5.  Cl.  and  Jr  arc  the  flux 
densities,  the  solid  jngles  into  which  the  flux  is  emitted,  and  the 
bandwidths  in  km  s " 1  for  OH  emission  and  35  jim  absorption 
respectively,  and  .-fj,  is  the  interstellar  attcntuation.  In  the 
simplest  case,  the  last  three  factors  are  1  and  the  predicted 
efficiency  for  a  saturated  maser.  given  by  the  flux  density  ratio. 
isO  25  (F.liizur.  Goldreich.  and  Seoville  1976).  — 

For  the  present  stars  Ssi  was  estimated  from  a  linear 
extrapolation  between  12  and  20  jim.  The  derived  ratio 
Sou  S,,  is  plotted  in  Figure  S  against  /.F,n  F,*.  Given  the 
uncertainties,  the  flux-density  ratio,  reflecting  the  efficiency,  is 
generally  close  to  the  expected  value  of  0.25.  Nevertheless, 
there  is  some  indication  of  a  trend  for  the  flux-densitv  ratio  to 
increase  from  0.15  to  0.6  with  increasing  circ  .  "i stellar 
reddening,  suggesting  a  corresponding  increase  in  the  -.creed 
pump  efficiency.  This  can  be  reconciled  with  the  expected  value 
c  =  0.25.  if  the  bandwidth  ratio  decreases  with 

increasing  thickness  of  the  circumstellar  shell.  When  the  shell 
becomes  thicker,  the  radius  at  which  the  majority  of  35  pm 


Fin  * — The  ratio  of  OH  to  J5  *im  flu  v  as  »  function  of  the  fraction  of  the 
infrared  energy  emitted  al  CO  tan. 


(4) 


"here  and  Fm„  mm  are  respectively  the  OH  and  bolo- 

metric  fluxes  at  the  maximum  and  minimum  of  the  light  curve. 
Thus,  the  relative  amount  of  35  pm  flux,  which  reflects  the 
temperature  of  the  outer  dust  shell,  should  change  as  a  function 
of  phase  in  the  light  curve.  Since  the  depth  of  the  silicate  feature 
indicates  a  temperature  gradient  between  the  regions  where  the 
10  jim  continuum  and  silicate  features  arise  (Kwan  and  Seo- 
vflle  19761.  r„,  should  also  vary  over  the  light  curve.  From 
equations  t2)  and  14)  it  can  be  shown  that 


(5) 


Thus,  combining  equations  (5)  and  13).  the  change  in  the  silicate 
optical  depth  between  minimum  and  maximum  of  the  light 
curve  is  given  by 


where  •«,,».  ms,  3r,=  'he  silicate  optical  depths  at  minimum  and 
maximum  light.  Because  of  the  large  width  of  the  10  jim  filter 
used  at  WIRO.  the  precise  value  of  A:  is  uncertain,  but  effec¬ 
tively  the  relation  implies  that  the  silicate  optical  depth  is 
smaller  at  maximum  light  than  at  minimum  and  should 
increase  with  decreasing  L.  This  is  in  good  qualitative  agree¬ 
ment  with  our  observations  of  those  stars  in  the  galactic-center 
sample  for  which  infrared  fluxes  were  measured  at  two  epochs 
(cf  Fig.  2  and  Table  6i  and  with  the  measurements  of  the  bright 
OH  IR  starOH  26.5  0.6 1  Forrest  etal.  197S). 


c)  Circumstellar  Shell  Si:e  unil  OH  Luminosity 
The  correlation  found  between  LOH  and  the  size  of  the  eir- 
cumstelljr  shell  by  BH  can  now  be  understood.  As  the  mass- 
loss  rate  increases,  the  circumstellar  shell  density  rises  so  that 
the  OH  maser  can  operate  at  larger  radii  from  the  star:  the 
minimum  density  required  for  a  saturated  maser.  nH  =  lO*  (cf. 
Elitzur.  Goldreich.  and  Scowlle  1976)  is  found  over  a  greater 
surface  area.  Since  the  stimulated -emission  rate  for  a  saturated 
maser  ts  determined  by  the  number  of  pump  photons.  Lu„  will 
not  be  enhanced  merely  by  enlarging  the  maser  volume. 
Nevertheless,  with  rising  density  the  circumstellar  dust 
becomes  cooler  and  more  stellar  emission  is  converted  to  35 
jim  flux.  Thus,  with  increasing  mass-loss  rate,  the  effects  of  a 
higher  number  of  pump  photons  jnd  of  available  OH  molecu¬ 
les  combine  to  produce  enhanced  maser  emission  and  keep  the 
maser  saturated. 


V.  CONCLUSION'S 

In  this  paper  wc  have  presented  simultaneous  OH  and  infra¬ 
red  measurements  of  j  homogeneous  sample  of  OH  IR  stars  at 
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the  tangential  paint  and  near  the  galactic  center:  their  dis¬ 
tances  are  therefore  well  known.  Very  accurate  radio  positions, 
derived  from  VLA  observations,  and  identifications  of  the 
infrared  counterparts  at  wavelengths  longward  of  8  p m  reduce 
the  possibility  of  confusion.  W  e  find  that  the  characteristics  of 
the  objects  in  each  sample  are  very  similar  in  most  respects. 
Although  there  is  some  tendency  for  the  luminosities  of  the 
galactic-center  stars  to  be  lower  than  those  at  the  tangential 
point,  the  paucity  of  objects  makes  it  difficult  to  ascertain 
whether  this  is  significant.  Future  studies,  perhaps  based  on  the 
results  of  the  IRAS  survey,  should  make  it  possible  to  verify 
w  hether  a  population  of  low -luminosity  masers  indeed  exists  in 
the  galactic  bulge. 

As  a  result  of  the  homogeneity  of  the  samples,  the  reliability 
of  the  infrared  identifications,  and  the  simultaneous  radio  and 
infrared  measurements,  it  has  been  possible  to  establish  a 
number  of  relations  between  various  stellar  and  circumstellar 
properties.  In  particular,  it  has  been  shown  quantitatively  that 
the  OH  luminosity  is  dominated  by  the  fraction  of  the  infrared 
energy  emitted  longward  of  20  nm  and  thus  by  the  stellar 
mass-loss  rate  However,  Lotl  is  to  some  extent  dependent  on 
the  boiometric  luminosity  so  that,  as  observed,  time  variations 
in  L„  are  reflected  in  LOH.  It  has  also  been  demonstrated  that 
the  apparent  optical  depth  of  the  silicate  feature  t„,  is  deter¬ 
mined  by  the  stellar  mass-loss  rate,  confirming  that  this  feature 
arises  in  the  circumstellar  shell.  Moreover,  using  the  depen¬ 
dence  of  both  LoH  and  r„,  on  the  stellar  mass-loss  rate,  and  the  - 
fact  that  the  amplitude  of  variability  of  the  OH  maser  emission 
is  half  that  of  the  stellar  emission,  the  observed  variations  in  r„, 
over  the  light  curve  of  a  given  object  can  be  explained  in  a 
qualitative  manner. 

Our  new  data  have  permuted  an  investigation  of  the  maser 
pump  efficiency.  Within  the  uncertainties,  this  is  found  to  be 
consistent  with  the  canonical  value  of  0.25.  However,  some 


tendency  for  the  inferred  efficiency  to  increase  with  increasing 
circumstellar  reddening  is  discernible  in  the  data. 

The  results  discussed  above  were  derived  from  a  sample  of 
stars  covering  a  wide  range  of  masses  and  luminosities. 
Although  these  objects  represent  the  high-luminosity  tail  of  the 
OH  luminosity  distribution,  it  seems  justifiable  to  extend  our 
conclusions  to  the  majority  of  OH  IR  stars.  In  particular,  it  is 
clear  that  sources  whose  stellar  properties  are  the  same  can 
have  a  variety  of  circumstellar  attributes.  In  effect,  for  a  given 
L,  a  large  range  of  values  of  LOH.  which  reflect  differing  mass- 
loss  rates,  is  possible.  Thus  these  stars  can  be  interpreted  as 
representing  an  evolutionary  sequence  of  increasing  mass-loss 
rate. 

The  measurements  described  here  were  made  at  six  observa¬ 
tories  in  the  course  of  seven  observing  runs.  As  a  result,  thanks 
are  due  to  many  people  for  their  contributions  to  the  success  of 
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APPENDIX  A 

DERIVED  PHYSICAL  PARAMETERS  FOR  OH  IR  STARS 

The  distance  to  the  OH  IR  stars  studied  here.  D,  is  based  on  a  galactic  distance  scale.  center-to-Sun.  of  R0  =  10  kpc.  For  the 
^tangential  sample  it  is  derived  kinematically  from  the  stellar  radial  velocity  r  =  (HV  +  LV)  2  Ikm  s ~ 1 )  iwhere  HV  and  LV  are  the 
velocities  with  respect  to  the  local  standard  of  rest  of  the  high-  and  low-velocity  components  of  the  OH  profile!,  using  the  galactic 
rotation  curve  of  Burton  and  Gordon  (-1978).  For  stars  with  radial  velocities  larger  than  the  maximum  velocity  allowed  by  circular 
rotation  rm„.  D  is  presumed  equal  to  the  distance  to  the  tangential  point.  Stars  with  r  <  r*,,  probably  lie  between  the  listed  near- 
and  far-kinematic  distances.  For  the  galactic-center  sample,  stars  within  a  degree  of  the  center  have  high  radial  velocities  1 1  r|  >  100 
km  s" ')  and  are  assumed  to  be  part  of  the  bulge  population  at  a  distance  of  10  kpctHabing  et  at.  1983).  At  I  >  1°,  where  the  density 
of  the  bulge  drops  off  significantly,  distances  are  more  uncertain.  Objects  with  small  radial  velocities  ( 1 1  |  S  10  km  s ~  *)are  probably 
foreground  stars  and  are  assumed  to  be  at  the  location  of  the  molecular  ring  at  5  kpc  distance,  w  here  the  density  distribution  of  the 
OH  IR  stars  peaks  icf.  BH.MW’l.  The  remainder  of  the  sample,  with  high  negative  or  positive  radial  velocities,  is  assumed  to  be  at  a 
distance  of  10  ±  10/  kpc.  where  /  is  the  source  galactic  longitude  in  radians. 

The  expansion  velocity  of  the  circumstellar  shell  is  r,  «  (HV  -  LV)  2  (km  s  *'). 

The  OH  luminosity  Lotl  at  1612  MHe  is  LOH  =  D4(Slv  w 1  (Jy  kpe:),  where  D  is  in  kpc  (cf.  BHMW'V  The  3  o  sensitivity  limit  of 
the  Dwingeloo  and  Onsala  survey,  from  which  the  tangential  sample  was  drawn,  is  '■0.6  Jy.  At  a  mean  distance  of  8  9  kpc  this 
corresponds  to  a  minimum  detectable  OH  luminosity  of  50  Jy  kpc:  and  to  a  completeness  limit  of  100  Jy  kpcJ.  For  the  galactic- 
center  sample,  similar  values  apply. 

Regarding  stellar  luminosity,  most  energy  distributions  in  Figures  1  and  2  rise  between  3  and  8.7  jim  and  many  appear  to  be 
_.^oader  than  expected  for  a  single-temperature  Planck  curve,  indicating  the  presence  of  circumstellar  dust  at  different  temperatures. 
The  stellar  radiation,  which  peaks  around  1-2  #im.  is  absorbed  by  the  dense  dust  shell  and  reemitted  at  longer  infrared  wavelengths. 
The  integrated  infrared  flux  Fm  determined  from  the  observed  energy  distribution,  and  including  a  correction  factor  of  25 "«  to 
account  for  flux  emitted  beyond  25  fim  (cf.  Werner  et  al.  1980),  is  therefore  a  good  measure  of  the  boiometric  luminosity  of  the  star. 


No.  2.  1985 


OH  1R  STARS 


639 


The  stellar  luminosity  is  then  L ,  =  2.9 D!F,t  (LQ).  where  D  is  in  kpc  and  F„  is  in  10  u  W  m‘*.  If  the  OH  sources  without  infrared 
detections  (Table  5|  have  energy  distributions  similar  to  those  detected,  their  10  pm  upper  limits  correspond  to  L  s  3.5  x  103  L 
for  a  mean  distance  of  8.9  kpc.  *  ° 

The  color  temperature  Tc  for  each  source  was  determined  by  fitting  a  Planck  curve  to  the  observed  infrared  energy  distribution 
between  4  8  and  12.5  j/m. 

As  for  silicate  absorption,  the  optical  depth  of  the  silicate  feature  t„,  seen  at  9.7  pm  in  Figures  1  and  2  was  derived  from  the  ratio  of 
the  measured  flux  at  the  deepest  point  of  the  feature  /.F„,  to  the  continuum  flux  at  the  same  wavelength  kFt„,.  interpolated  from  the 
8.7  and  the  12.6  pm  flux  measurements,  so  that  rtil  =  — In  (kFvVkFM J. 


APPENDIX  B 

THE  EFFECTS  OF  INTERSTELLAR  EXTINCTION 

Stellar  parameters,  such  as  the  bolometric  luminosity  and  the  color  temperature,  which  are  derived  from  the  infrared  observa¬ 
tions.  will  be  influenced  by  interstellar  extinction.  This  can  range  from  7  to  32  mag  in  the  visual  for  these  distant  stars  (de  Jong  1983). 
Using  the  extinction  law  in  the  infrared  from  Beckliner  of.  11978),  we  calculate  that  for  an  average  visual  extinction  of  20-30  mag  the 
observed  luminosity  of  a  400  K  blackbody  will  be  underestimated  by  a  factor  of  ~  1.6.  Since  both  distance  uncertainties  and 
variations  in  the  extinction  with  respect  to  the  mean  could  cause  uncertainties  in  the  derived  bolometric  luminosities  of  a  factor  of 
2-3  for  individual  objects,  no  corrections  for  interstellar  extinction  were  made.  However,  the  influence  of  25  mag  of  extinction  on  a 
blackbody  of  400  K  is  show  n  in  Figure  3  by  an  arrow.  None  of  the  trends  discussed  in  the  body  of  the  text  appear  to  be  significantly 
influenced  by  interstellar  extinction. 
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INFRARED  OBSERVATIONS  OF  THE  EARLY  UNIVERSE 
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Abstract:  The  possibilities  for  making  infrared 
observations  of  great  scientific  interest  from  above 
the  Earth' 8  atmosphere  are  discussed.  It  is  shown 
that  a  telescope  system  like  the  Large  Aperture 
Infrared  Telescope  System  (LAIRTS),  now  undergoing 
development  by  the  Air  Force  should  be  able  to  make 
certain  unique  observations  which  are  beyond  the 
capabilities  of  presently  existing  telescopes.  These 
observations  include  galaxies  in  the  process  of 
formation  in  the  early  Universe,  and  searches  for 
brown  dwarf  components  of  galactic  halos. _ 


INTRODUCTION 

The  latter  half  of  the  20th  century  has  seen  the 
introduction  of  a  large  number  of  new  techniques  in 
observational  astronomy,  a  revolution  which  continues 
today.  The  revolution  has  opened  wide  the  entire 
electromagnetic  spectrum  for  astronomical  observa¬ 
tions.  Visual  observations,  utilizing  photographic 
and  photoelectric  methods  developed  during  the  late 
19th  and  early  20th  centuries,  are  being  complemented 
by  ground-based  radio  and  Infrared  observations;  and 
satellite  observatories  operating  at  ultraviolet,  x- 
ray,  and  gamma-ray  frequencies  are  now  realities.  Our 
knowledge  of  the  universe  has  also  been  substantially 
enhanced  by  high  altitude  observations  at  non-optical 
frequencies  from  aircraft,  balloons,  and  sub-orbital 
rocket  flights.  In  1983  the  Infrared  Astronomy 
Satellite  (IRAS)  completed  a  high  sensitivity  survey 
over  9555  of  the  entire  sky,  demonstrating  the  feasi¬ 
bility  of  making  long-term  infrared  observations  from 
space . 

A  new  instrument  conceived  and  being  developed  by  a 
group  at  the  Air  Force  Geophysical  Laboratory  (AFGL) 
in  Massachusetts  is  the  Large  Aperture  Infrared 
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Telescope  System  (LAIRTS).  As  envisioned,  it  will 
provide  man  with  the  next  significant  step  forward  in 
Infrared  observational  capability.  LAIRTS  is  to  be 
launched  above  Barth's  atmosphere,  with  the  entire 
telescope  cooled  to  S25K.  Its  major  advantage  over 
ground-based  Infrared  telescopes  will  be  its 
sensitivity,  which  depends  on  its  cooled  optics,  and 
its  location  above  the  absorbing  and  emitting 
atmosphere.  For  observations  of  extended  sources,  and 
for  mapping  of  pre-selected  regions  of  space,  LAIRTS 
will  be  fitted  with  an  Infrared  imaging  camera,  making 
possible  high  sensitivity  measurements  at  high  spatial 
and  medium  spectral  resolution.  Certain  observations, 
to  be  discussed  herein,  cannot  be  made  at  all  by 
presently  existing  telescope  systems,  nor  by  others 
expected  to  become  operational  in  the  near  future, such 
as  the  Hubble  Space  Telescope.  LAIRTS  will  thus  be  an 
Important,  and  presently  unique,  tool  for  making  these 
types  of  observations. 


LAIRTS  IMAGING  CAMERA 
Sensor  Description 

The  LAIRTS  telescope  design  is  a  four  mirror  off- 
axis  folded  Schmidt,  with  no  obscurations  or  spiders. 
The  primary  mirror  is  a  32"  x  22"  ellipse,  with  3670 
cm2  collecting  area.  The  optical  system  is  diffrac¬ 
tion  limited  for  all  X  >  S  iim. 

A  dichroic  beamsplitter  reflects  70*  or  more  of  |i  < 
30  iim  radiation  toward  the  three-mirror  camera  system, 
which  then  re-images  the  telescope  field  onto  a 
detector  array.  The  array  consists  of  four  64  x  64 
Si: As  direct  readout  mosaics,  arranged  in  a  square. 
The  field  of  view  (F 0V)  for  each  pixel  is  an  area  on 
the  sky  11  arcsec  x  11  arcsec  square,  with  12  arcsec 
center-to-center  distance  between  adjacent  pixels. 
Adjacent  64  x  64  mosaics  are  separated  by  48  arcsec, 
so  that  the  total  array  covers  a  field  26.4  arc- 
minutes  square. 

Sensitivity 

The  sensitivity  of  the  LAIRTS  Imaging  Camera  is 
expected  to  be  background  limited  by  zodiacal  emission 
(ZB).  A  model  for  the  ZB,  adopted  from  the  SIRTF 
Interim  Report  (1978),  was  used  to  calculate  flux 
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density  from  the  ZE  which  the  LA1RTS  camera  will 
detect  at  selected  wavelengths.  The  results  are  shown 
In  Table  1.  These  estimates  agree  reasonably  well 
with  ZE  brightness  observations  made  by  IRAS  at  12  and 
2S  pm  In  the  direction  of  the  ecliptic  pole  (Hauser  et 
al.,  1984). 


Table  1.  Zodiacal  Emission 

X ( urn)  Flux  Density  (mJv) 

5  7.6 

10  57 

20  38 

The  noise  equivalent  power  (NEP)  for  each  element  of 
an  array  detector  is  given  by  Lamb  et  al.  (1984)  as 

NEP  =  E.  -  [2(n*+e) ]1/2WHz_1/2  (1) 

A  r> 

where  Ex  =  photon  energy 

r>  =  quantum  efficiency 

9  =  photon  flux  (photons/ sec) 

$  =  dark  current  (electrons/sec) 

Since  the  camera  is  expected  to  be  background  limited 
by  the  photon  flux  from  the  ZE,  its  NEP  may  then  be 
calculated  from 


NEP  *  Ex 


1/2 

WHz 


-1/2 


(2) 


where  r»  =  0.3. 

The  lo  sensitivity  for  an  integration  time  ^  =  l 

second  is  then 


lo 


sensitivity  — 


NEP 


W  cm-2 


(3) 
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where  Ac  =  collector  area  =  3570  cm2 

£  =  system  transmission  due  to  optical  losses 

»  0.49  -  (0.98)7(0.80) (0.70) 

[7  mirrors] [filter] [dichroic] 

Using  equations  (2)  and  (3)  the  noise  equivalent  flux 
densities  (NEED'S)  of  the  camera  for  detection  of  both 
point  and  extended  sources  were  calculated  for 
selected  wavelengths  and  spectral  bandwidths.  The 
result  of  these  calculations  are  presented  in  Tables  2 
and  3 . 

Table  2.  NEFD  of  Imaging  Camera  for  point  source 
(t,  *  1  second) 

X ( um)  AX  °  0. IX  AX  =  0.5X 

5  440  pJy  200  |iJy 

10  1200  |lJy  540  pJy 

20  980  pJy  440  pJy 


Table  3.  NEFD  of  Imaging  Camera  for  extended  source 


( t  ^  =  1  second ) 


X  ( pm) 

AX  =  0.1X 

AX  «  0.5X 

5 

3.6 

nJy/sq-arcsec 

1.7 

p Jy/sq-arcsec 

10 

9.9 

p Jy/sq-arcsec 

4.5 

p Jy/sq-arcsec 

20 

8.1 

M Jy/sq-arcsec 

3.6 

p Jy/sq-arcsec 

UNIQUE  OBSERVATIONS 
Objects  of  High  Redshlft 

All  astrophysical  sources  at  great  distances 
exhibit  doppler-shif ted  spectra,  indicating  velocities 
of  recession  proportional  to  their  distance.  The 
observed  wavelengths  of  radiation  are  .shifted  toward 
the  red  from  the  emitted  (or  rest)  wavelengths.  The 
corresponding  cosmological  redshlft,  Z,  is  defined  in 
terms  of  the  observed  wavelength,  X',  and  the  emitted 
wavelength,  X: 
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Z  - 


OX 

x~ 


The  relationship  between  recession  velocity,  v,  and 
redshift,  Z,  is: 


Z 


r  i+v/c  1/>2 

l l-V/C 


It  follows  that  observed  energy  distributions  from 
distant  sources  are  also  redshlfted.  Contributions  to 
the  energy  distribution  redshift  arise  from  3  effects: 

( 1 )  The  photon  energy  decreases 


E 


v 


V 

1+Z 


(2)  A  redshifted  frequency  interval  becomes 
smaller 


ov ' 


ov 

1+Z 


(3)  The  rate  of  photon  emission  decreases  due  to 
time  dilation 


A 


t 


The  net  result  of  these  effects  on  the  energy 
distribution  is  to  decrease  the  observed  flux  density 
according  to  the  relation 

Fy 1  =  KF„W  M-2Hz-l 

where 


r 


2(  1+Z) 
(1+Z)2+1 


K 


1 


1/2 
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Values  of  K  and  recession  velocities,  V,  for  selected 
redshlfts,  Z,  are  given  in  Table  4. 

Table  4.  Values  of  K  and  V  for  Selected  Values  of  Z 

Z  V  K 


0.3 

0.257c 

0.967 

1.0 

0.600c 

0.800 

2.0 

o.aooc 

0.600 

5.0 

0.946c 

0.324 

10 

0.984c 

0.180 

20 

0.9955c 

0.0950 

100 

0.9998c 

0.0198 

Giant  elliptical  galaxies  are  among  the  brightest 
visible  objects  known,  making  them  among  the  easiest 
objects  of  high  redshift  to  detect.  A  model  galaxy 
was  constructed  based  on  the  energy  distribution 
discussed  by  Grasdalen  (1980).  The  photometry, 
colors,  and  absolute  K  magnitudes  published  by 
Lebofsky  (1980)  were  used  to  set  an  absolute  scale  for 
Grasdalen* s  energy  distribution.  The  absolute  visual 
magnitude  of  the  model  is  *>  -23.1.  The  apparent 
visual  magnitude  is  then  m.,  «  0.0  at  a  distance  of 
0.41  Mpc,  too  close  for  an  object  to  exhibit  a 
measurable  cosmological  redshift.  Spectral  energy 
flux  densities  were  then  calculated  for  the  model 
galaxy  at  selected  values  of  the  redshift,  Z, 
including  effects  of  the  inverse-square  law  for 
radiant  emission  by  point  sources.  At  distances 
corresponding  to  Z  >  1.0,  all  galaxies  will  appear  as 
point  sources  in  the  field  of  the  LAIRTS  imaging 
camera.  No  model  for  galaxy  evolution  was  adopted, 
and  no  contribution  due  to  a  possible  non-zero  value 
for  the  cosmological  deceleration  parameter  (qQ)  was 
included.  The  resulting  expected  flux  densities  for 
giant  elliptical  galaxies  at  various  redshifts  are 
given  in  Table  5.  An  examination  of  the  table  shows 
that  the  wavelength  of  peak  emission  is  shifted 
continually  red-ward  for  Increasing  Z. 

The  model  predicts  that  a  normal  giant  elliptical 
galaxy  observed  at  Z  e  5  will  exhibit  a  peak  flux 
density  of  -55pJy  at  10  pm.  With  an  NEPD  of  540  pjy 
for  1  second  integration  (and  aX  *  Sum) ,  a  40-minute 
integration  time  would  be  required  to  obtain  a 
signal/noise  ratio  of  5.  For  Z  =  10  the  peak  emission 
is  near  20  pm,  with  a  value  of  B28  pjy.  In  this  case. 


Table  S.  Expected  Flux  Densities  for  Giant  Elliptical  Galaxies  at  Various  Redshlfts  (Z) 
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a  detection  with  S/N  =  6  (and  oX  =6  pm)  would  need 
146  minutes  of  integration  time. 

Primordial  Galaxies 


Of  particular  interest  are  galaxies  with  redshifts 
of  Z  *  5-10.  According  to  one  view  of  the  early 
universe,  proposed  by  Zel'dovlch  (1973),  galaxies  were 
formed  by  condensation  from  initial  mass  perturbations 
the  size  of  galaxy  clusters  or  superclusters. 
Galaxies  thus  formed  are  expected  to  collapse  at 
redshifts  -5  or  less,  and  no  galaxies  with  Z  substan¬ 
tially  greater  than  5  would  be  expected  to  be 
observed.  A  distinctly  different  scenario,  due  to 
Press  and  Schechter  (1974),  depicts  galaxies  being 
built  up  from  smaller,  globular  cluster  size  units. 
In  this  case  galaxies  are  expected  to  form  much 
earlier,  at  Z's  of  20  or  greater. 


According  to  at  least  some  galaxy  formation  models 
(c.f.  Kaufman  1976),  galaxy  birth  should  be  accom¬ 
panied  by  a  burst  of  star  formation,  with  a  period  of 
peak  brightness  lasting  -2  x  107  years,  during  which 
the  galaxy  luminosity  should  be  dominated  by  hot  0  and 
B  star  emission.  This  should  make  the  spectral  energy 
distribution  bluer  than  our  galaxy  model  of  Table  5, 
and  about  a  factor  of  100  brighter.  The  predicted 
integration  times  to  detect  bright  starburst 
primordial  galaxies  at  various  redshifts,  with  ax  = 
0.1X,  and  S/N  >  10  are  given  in  Table  6. 


Table  6. 


t,  (sec)  to  Detect  Primordial  Galaxies; 
S/N  >  10;  OX  =  0.1X 


X  ( urn) 


z  =  5  z  -  10  z  =  20 


5 

10 

20 


3  196 

5  68  3700 

14  147 


Massive  Galactic  Halos 


Recent  dynamical  evidence  (Einasto  et  al . ,  1974; 
Bok  1983)  suggests  that  our  Galaxy  contains  -5  x 
101:1Me  of  Invisible  matter  in  an  extensive  halo 
(radius  >  60  kpc)  surrounding  the  readily  observed 
visible  disc.  Similar  evidence  for  the  existence  of 
massive  halos  has  been  obtained  for  a  number  of 
external  spiral  galaxies  (Roberts  1975;  Rubin  et  al . , 
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1978).  A  list  of  these  galaxies  Is  presented  In  Table 
7. 


Table  7.  Spiral  Galaxies  with  Flat  Rotational  Curves 


NGC  4378 
NGC  4594 
NGC  7217 
NGC  2590 


NGC  1620 
NGC  3145 
NGC  801 
NGC  7541 


NGC  7664 
NGC  2998 
NGC  3672 


There  is  at  present  no  observational  evidence  to  give 
us  a  clue  as  to  the  form  or  physical  properties  of 
these  massive  invisible  halos  (Boh  1983).  Suppose, 
however,  that  such  a  halo  was  completely  made  of  brown 
dwarf  stars,  those  theoretical  objects  midway  in  size 
between  Jupiter  and  the  sun,  not  quite  massive  enough 
to  ever  have  initiated  nuclear  reactions  in  their 
cores.  These  stars  would  be  cool  and  would  emit 
energy  at  infrared  frequencies.  A  possible  model  for 
a  brown  dwarf  galactic  halo  could  have  the  following 
characteristics : 


Halo  Shape: 

Halo  Radius: 

Halo  Mass: 

Brown  Dwarf  Radius: 

Brown  Dwarf  Mass: 

Total  No.  of  Brown  Dwarf  Stars.- 
Distance  to  Spiral  Galaxy: 


spherical 
60  kpc 
5  x  1011mq 
0.1  Rq 
0.01  Mo 

5  x  1013 
30  Mpc 


The  diameter  of  a  spiral  galaxy  the  same  size  as 
the  Milky  Way  (radius  -  15  kpc)  at  a  distance  of  30 
Mpc  would  subtend  an  angle  of  -3.4  arc-minutes.  A  60 
kpc  radius  halo  would  then  lie  inside  20  arc-minutes, 
and  would  be  included  in  the  26  x  26  arc-min  field  of 
the  LAIRTS  camera.  Based  on  the  model  parameters 
listed  above,  the  predicted  flux  densities  which  would 
need  to  be  measured  by  the  LAIRTS  Imaging  Camera  are 
given  in  Table  8  for  brown  dwarfs  of  various  tempera¬ 
tures.  For  convenience,  the  calculated  system  NEFD  is 
also  given.  Note  that  brown  dwarf  halos  down  to 
temperatures  of  “300  K  can  be  observed  at  10  Mm  with 
S/N  -10  in  less  than  5  minutes,  provided  the  model  is 
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basically  accurate.  Such  observations  will  be 
important.  Even  null  results  will  be  scientifically 
interesting,  since  they  will  set  certain  limitations 
on  the  form  and  characteristics  of  the  matter  which 
makes  up  the  massive  invisible  halos. 

Table  8.  Expected  Flux  Densities  from 
Brown  Dwarf  Halos 


( ti  Jy/square-arcsecond) 


(»M) 

T=150K 

T=300K 

T=600K 

T=1000K 

NEFD* 

5 

1 .OxlO-4 

1 . 5 

189 

1350 

3.6 

10 

0.48 

5 . 9 

71 

220 

9.9 

20 

.  18 

2 . 2 

9.5 

21 

8.1 

SUMMARY 

The  potential  for  making  infrared  observations  of 
great  scientific  interest  exists.  Observations  of 
galaxy  formation  in  the  early  universe  should  be 
possible  using  LAIRTS,  which  is  now  under  development 
by  the  Air  Force  Geophysical  Laboratory.  Observa¬ 
tions  of  brown  dwarf  components  of  galactic  halos  is 
also  shown  to  be  a  possibility  with  LAIRTS.  The 
implementation  of  both  of  these  observing  programs 
will  significantly  expand  our  knowledge  of  the 
universe  in  which  we  live. 
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ABSTRACT 

Infrared  photometry  of  Nova  Vulpeculae  1984  number  2  (NV2)  from  2.3  to  19.5  pm  during  1985  Mav  14  to 
1986  March  31  shows  that  silicate  grains  had  condensed  and  grown  in  the  nova  ejecta  by  1985  August  23.  240 
days  after  the  eruption.  A  relative  overabundance  of  oxygen  in  the  nova  shell  seems  indicated.  Forbidden  12.8 
p m  [Ne  n]  emission  was  a  factor  of  *  41  above  the  continuum  at  a  spectral  resolution  of  A/AX  =  67  on  dav 
240:  the  line  persisted  through  day  461.  The  anomalous  chemical  composition  of  NV2's  ejecta  supports  recent 
suggestions  that  ONeMg  white  dwarfs  evolved  from  8-12  M  ,  progenitor  stars  are  accreting  matter  m  binary 
systems. 

Subject  headings:  infrared:  sources  —  stars:  individual  —  stars:  novae 


I.  INTRODUCTION 

We  report  the  discovery  of  the  10  and  20  jim  silicate 
emission  features  in  Nova  Vulpeculae  1984  number  2  (NV2) 
and  suggest  that  silicate  grains  have  formed  in  the  ejecta.  The 
observations  provide  an  estimate  of  the  total  mass  of  silicate 
grains  condensed  in  the  shell. 

The  observations  reported  here,  together  with  our  earlier 
discovery  in  NV2  of  strong  [Ne  aj  emission  at  12.8  pm 
(Gehrz,  Grasdalen,  and  Hackwell  1985:  hereafter  Paper  I), 
confirm  the  presence  of  oxygen-neon-magnesium  (ONeMg) 
white  dwarfs  accreting  matter  in  binarv  systems.  Such  novae 
may  be  capable  of  contributing  significant  quantities  of  the 
interesting  radioactive  isotopes  22Na  and  :6A1  to  the  interstel¬ 
lar  medium  (see,  for  example,  Hillebrandt  and  Thielemann 
1982  and  Truran  1984).  Meteoritic  abundance  anomalies  sug¬ 
gest  that  these  extinct  isotopes,  the  by-products  of  a  transient 
nucleosvntheiic  event  such  as  a  nova  eruption,  were  injected 
into  the  primitive  solar  system  during  the  early  stages  of  its 
formation  (Truran  1984). 

II.  OBSERVATIONS 

Photometric  and  spectrophotometric  observations  (Table  1 
and  Figs.  1-3)  were  obtained  with  the  Wyoming  bolometer 
(Gehrz.  Hackwell.  and  Jones  1974)  and  the  A/AA  -  67  As:Si 
spectrophotometer  (Gehrz  et  al.  1984)  on  the  234  cm  Wyo¬ 
ming  Infrared  Telescope.  Several  apertures  (5"  or  6"8)  and 
beam  separations  (10"- 18")  were  used,  the  measurements  are 
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!  Wyoming  Infrared  Observatory.  Department  of  Physics  and  Astron¬ 
omy.  University  of  Wyoming. 

'Space  Sciences  Laboratory.  The  Aerospace  Corporation.  Los  Angeles. 
California. 

4  Department  of  Physical  Sciences,  Eastern  Montana  College.  Billings. 
Montana. 


independent  of  these  parameters  because  NV2's  shell  is  spa¬ 
tially  unresolved.  Most  observations  reported  here  were  made 
with  remote  telephone  links  enabling  observers  to  control  the 
Wyoming  telescope  and  its  instrumentation  from  laboratories 
in  Laramie.  Wyoming,  and  Minneapolis.  Minnesota. 

III.  INFRARED  SILICATE  EMISSION 

Broad-band  photometry  and  narrow-band  spectrophotome¬ 
try  for  1985  August  23.4  show  the  presence  of  the  10  pm 
silicate  emission  feature  (Fig.  2).  The  20  pm  silicate  emission 
feature  was  present  on  August  23.4,  but  the  19.5  pm  flux  level 
on  May  15.4  was  consistent  with  the  free-free  continuum 
extrapolated  from  data  at  shorter  wavelengths  (Fig.  1). 

Both  silicate  features  should  have  been  detectable  over  the 
free-free  continuum  on  May  15.4  had  they  been  present  at  the 
level  measured  in  August  (see  Figs.  1-3).  The  data  imply  that 
the  silicate  emission  features  had  increased  in  intensity  by  at 
least  a  factor  of  2  between  May  15.4  and  August  23.4  due  to 
grain  condensation  and  growth. 

An  alternative  source  of  the  silicate  emission  could  be  the 
illumination  after  day  140  of  a  dust  shell  remnant  of  the 
mass-loss  phase  of  the  nova  progenitor  (see  Bode  and  Evans 
1983).  This  seems  unlikely.  First,  the  material  would  have  to 
lie  >  3.6  X  1017  cm  from  the  binary  system  and  would  sub¬ 
tend  an  angular  diameter  S  16".  NV2  is  unresolved  by  a  5" 
beam  at  10  pm.  Second,  the  temperature  for  grains  directly 
illuminated  by  starlight  at  this  distance  from  a  source  with 
NV2’s  outburst  luminosity  of  *  10'  L~,  (Paper  1.  erratum) 
will  be  £  40  K.  The  observed  grain  temperature  is  200-400 
K  on  day  240  assuming  that  the  20  pm/ 10  pm  opacity  ratio 
for  silicate  grains  lies  between  0.3  and  1.  We  cannot  entirely 
rule  out  the  possibility  that  the  progenitor's  ejecta  are  aniso- 
tropically  distributed  and  that  very  small  dust  grains  lying 
along  the  line  of  sight  of  our  beam  are  impulsively  heated  by 
hard  photons  (Desert.  Boulanger,  and  Shore  1985). 
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\F.  IN  WATTS  /cm 


TABLE  1 

Photometry  of  Nova  Vulpeculae  19X4  Number  2 
Magnitude 


A  Dav  140  Dav  171  Day  230  Dav  240  Day  294  Day  335  Day  3X2  Day  461 

<am>  (05/14/85)  (06/13/85)  (08/13/85)  (08/23/85)  (10/16/85)  (11/26/85)  (01/12/86)  (03/31/86) 


2.3 .  A- 7.29  ±  0.04  +7,42  +  003  +7.86  +  0.08  +  811  ±  0.03  +  8.30  +  0.05  ...  +8.62  ±  0.12  +  8  90  ±0  21 

3.6  .  +  6  09  +  0.02  +  6.22  ±  0.01  +  7.04  ±  014  +  6  87  +  003  -  6.63  +  0.05  ...  a-  6.93  ±010  -  6.59  -  0  06 

4  9 .  -  5.71  ±  0.08  +  5  61  ±  0.07  +  6.00  ±  0.24  +  6.34  ±  0.09  . .  *  6  91  ±  0.28  +  7.03  ±  0.15 

8.7  .  +4.43  ±  0.20  +4  03  ±  0  09  +  4  49  ±  0.33  -4.70  ±0.16  +4.44  ±  0.27  ...  +4.61  +  0.23  -4  34  +  0  14 

10.0a .  ...  +3.07  +  0.09  ...  +2.99  ±  0.05  -3.15  +  0.16  -  3.28  ±  0.02  +3  38  +  0  14  -3.28  +  008 

V"  .  +2.70  ±  0.07  +3.04  ±  0  05  +3.10  ±  008  -3.20  -  006  +3.30  ±  0.04  +2.67  ±  0.08  -  3.39  +  0.08  -  3.30  -  0  12 

11  4 .  +3.20  -  0.23  +  3  04  ±  0.10  +  2.76  -  0.28  -  2.81  +  0.03  -  2.91  +  0.08  -  2  02  +  0.32  -  3  44  -  0  20  -  3  19  +  0  16 

12.6 .  -  1.69  +  0.07  -2  05  ±  0  04  +2.25  ±014  -2.33±  0.08  +  2.38  +  0.06  -  1.95  ±  0.11  -  2.20  ±  0 10  -->95^018 

Qb .  +2.55  ±  0.29  +  1  97  ±  0.25  +  1.97  -  0.25  2- 1.75(3  a)  +1.83  ±  0  33  ...  +2.03  ±  0  20  -2  42  ±014 


WAVELENGTH  IN  MICRONS 

Fig  1 


WAVELENGTH  IN  MICRONS 
Flo.  2 


Fig  1  -Broad-band  infrared  energy  distributions  of  NV2  before  ( dasned  line )  and  after  ( solid  line)  the  formation  of  circumstellar  silicate  grams 
Narrow-band  data  ire  plotted  for  the  12  8  am  |Ne  it)  emission  peak  Free-free  lines  indicate  tF,  -  constant  and  have  been  normalized  to  photometry  at 
A  2  5  am  Hydrogen  (5  —  4)  was  present  in  the  3  pm  band  in  1985  May  (Gchra.  Grasadalen.  and  Haclcwell  1985).  Error  bars  are  smaller  than  the  plotting 
symbols  unless  indicated  The  neon  line  contaminates  the  114  am  passband  Both  the  10  and  20  am  silicate  features  are  present  in  1985  Augusi 

Fig.  2  —Infrared  spectrum  of  NV2  from  8  to  13.2  am  showing  the  12  8  am  |Ne  n)  emission  line  and  the  10  am  silicate  emission  feature  Triangles 
denote  broad-band  photometry  Vega  denotes  the  aero  magnitude  calibration  F.rrors  as  indicated  in  legend  for  Fig  1 


NEON  NOVA  II. 
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Fig.  3.—  Comparison  of  the  10  ixm  spectra  in  1985  May  ( Unshed  line) 
and  August  ( solid  line)  showing  that  the  silicate  emission  feature  should 
have  been  detected  well  above  the  free-free  con.muum  in  Mav  if  it  had 
been  present  at  the  level  measured  in  August.  Error  bars  are  as  indicated 
in  the  legend  to  Fig.  1. 


IV.  ABUNDANCES  OF  OXYGEN,  NEON,  AND 
GRAIN  CONSTITUENTS 

The  formation  of  silicate  grains  suggests  that  the  ejecta  of 
NV2  are  oxygen  rich.  It  is  believed  (Hackwell  1971, 1972)  that 
amorphous  carbon  and  silicon  carbide  (SiO  grains  form  when 
C/O  >  1,  and  that  silicate  grains  condense  when  C/O  <  l. 
Thus,  the  formation  of  carbon-rich  (Gehrz  et  al.  1980a;  Gehrz 
et  at.  19806)  and  SiC-rich  (Gehrz  etal  1984)  shells  in  novae 
presumably  results  from  eruptions  on  carbon-oxygen  (CO) 
white  dwarfs  evolved  from  low-mass  stars. 

In  Paper  I,  we  reported  the  discovery  of  forbidden  12.8  pm 
(Ne  u]  emission  from  NV2  with  a  peak  intensity  of  *  60  Jy 
in  a  0.2  pm  bandpass  on  1985  May  15.4.  We  argued  that  neon 


was  overabundant  in  the  ejecta  of  NV2.  C.  H.  Townes  (private 
communication)  continued  the  presence  of  the  (Ne  nj  emis¬ 
sion  feature  m  1985  June  The  12.8  am  (Ne  ll|  emission 
feature  had  a  peak  intensity  on  August  23.4  of  18.5  Jy  m  a  0.2 
pm  bandpass.  The  decline  in  the  line  strength  over  the  past 
vear.  considering  the  effects  of  shell  expansion,  suggests  that 
excitation  conditions  have  not  changed  appreciably.  Our  con¬ 
clusion  (Paper  I)  that  the  amount  of  neon  in  the  form  of 
[Ne  n)  alone  must  be  at  least  solar  abundance  still  obtains. 
Recent  optical  and  ultraviolet  spectroscopy  of  NV2  bv  the 
HE  nova  team  (S.  <i  Starrfield.  private  communication) 
reveals  strong  lines  of  Ne  in  and  Ne  tv  ions  as  well.  The 
infrared  and  IUE  data  demonstrate  that  neon  is  enhanced  in 
the  ejecta  of  NV2.  Starrtieid.  Sparks,  and  Truran  (1986i  also 
report  strong  UV  and  optical  emission  from  Mg  II  Mag¬ 
nesium.  predicted  to  be  enhanced  in  the  ejecta  of  ON'eMg 
white  dwarfs,  is  a  primary  constituent  of  the  minerals  which 
produce  the  10  and  20  pm  emission  features. 

We  can  estimate  a  lower  limit  to  the  total  mass  of  the 
silicate  material  in  the  ejecta  of  NV2.  The  ratio  of  the  Dux 
observed  from  the  10  pm  emission  peak  to  the  flux  that  would 
be  emitted  by  an  optically  thick  dust  shell  reradiating  the 
outburst  luminosity  of  *  I05  L,  (Paper  I,  erratum)  gives  the 
mean  shell  silicate  optical  depth  (r10).  The  total  mass  Vf  of 
silicates  in  the  shell  is 

Mv  *  ir«;(T10)/Ks,, 

where  R,  is  the  shell  radius  and  k  ~  3  X  10s  cmr  g~ 1  is  the 
10  pm  silicate  opacity.  Assuming  an  expansion  velocity  of 
K,  *  1000  km  s  ’1  for  the  dust-forming  component  of  the 
ejecta  (Paper  I),  the  shell  should  have  attained  a  radius  of 
2.1  x  I015  cm  by  day  240.  and  an  optically  thick  sheil  re¬ 
radiating  10'  L  ,  on  that  day  would  have  had  a  temperature 
of  =  600  K..  The  10  pm  flux  predicted  for  such  a  shell 
compared  to  the  intensity  of  the  10  pm  silicate  feature  on 
day  240  gives  (t10)  *  4  X  10"3  and  *  10"*  M-.  We 
note  that  the  hydrogen  mass  associated  with  the  silicate  dust 
alone  is  *  250  Mv  *  3  x  10"'’  M ->  if  the  Si  abundance  is 
normal.  This  is  a  fair  fraction  of  the  mass  loss  predicted  for 
an  explosion  on  the  surface  of  an  ONeMg  white  dwarf 
(StaiTfield,  Sparks,  and  Truran  1986)  and.  of  course,  excludes 
the  mass  associated  with  the  high-velocity  ejecta  (see  Paper  I). 
Our  calculations  underestimate  the  value  of  (tio)  if  the  nova 
remnant  subsided  to  the  Eddington  limit  shortly  after  out¬ 
burst  (see  Starrfield.  Sparks,  and  Truran  1986)  and  if  the  nova 
remnant  luminosity  declined  from  the  Eddington  value  there¬ 
after.  Thus,  our  estimate  for  M„  represents  a  lower  limit. 

v.  DISCUSSION 

We  have  argued  that  the  silicate  emission  in  NV2  indicates 
that  grains  nucleated  and  grew  in  the  slow  ejecta  of  NV2. 
Models  of  runaway  nuclear  explosions  on  the  surface  of  an 
ONeMg  white  dwarf  show  that  about  half  the  ejected  mass 
travels  outward  at  velocities  as  low  as  500-1000  km  s 
which  is  consistent  with  the  rather  long  time  scale  on  which 
NV2  was  observed  to  condense  dust  grains.  Because  silicates 
have  a  very  low  etnissivity  in  the  near- infrared,  the  late 
appearance  of  the  shell  is  not  unexpected  and  is  generally 
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consistent  with  the  proposition  that  the  grains  nucleated  and 
grew  after  the  ejecta  in  which  the  grains  condensed  passed  the 
1000  K  radiation  field  point  in  the  expanding  shell  (Gallagher 
1977:  Gehrz  et  at.  1980a;  Gehrz  eta/.  198 Ob).  Starrfield, 
Sparks,  and  Truran  (1986)  have  suggested  that  the  anomalous 
emission  from  the  ejecta  of  NV2  can  be  explained  by  a 
runaway  nuclear  explosion  on  the  surface  of  an  ONeMg  white 
dwarf.  Law  and  Ritter  (1983)  suggested  that  some  nova 
systems  might  contain  ONeMg  white  dwarfs,  and  Nomoto 
(1983)  has  shown  theoretically  that  these  white  dwarfs  could 
be  the  remnants  of  8-12  Mc  main-sequence  progenitors. 

The  data  presented  in  Paper  1  and  by  Starrfield,  Sparks, 
and  Truran  (1986)  establish  that  two  of  the  crucial  ingredi¬ 
ents.  neon  and  magnesium,  are  enhanced  in  the  ejecta  of 
NV2.  Our  discovery  of  the  existence  of  silicate  grains  in  the 


ejecta  of  NV2  shows  that  there  is  an  enhancement  of  the  third 
crucial  ingredient  of  an  ONeMg  white  dwarf,  oxygen,  in  this 
nova. 


We  thank  J.  G.  Cohen,  J.  S.  Gallagher,  F.  C.  Gillett.  E.  P. 
Ney.  G.  W.  Preston,  S.  G.  Starrfield.  and  J.  W.  Truran  for 
stimulating  discussions.  C.  H.  Townes  communicated  his 
spectra  of  NV2  prior  to  publication.  L.  D.  Chisholm.  L.  R. 
Shaw,  C.  Jaworowski.  and  T.  Williams  assisted  with  the 
observations  and  maintained  the  Wyoming  systems.  R.  D.  G. 
is  supported  by  the  Institute  of  Technology  at  Minnesota. 
Astronomy  at  Wyoming  is  funded  by  the  National  Science 
Foundation,  the  USAF,  and  the  Department  of  Physics  and 
Astronomy. 
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A  SEARCH  FOR  COOL  COMPANIONS  OF  PLANETARY  NEBULA  NUCLEI 


A.  F.  Bentley 

Department  of  Physical  Sciences 
Eastern  Montana  College 
Billings,  Montana  59101 


At  present  only  a  small  number  of  planetary  nebulae  are  known  to  possess 
binary  nuclei.  Since  approximately  2/3  of  main  sequence  stars  are  mem¬ 
bers  of  binary  or  multiple  star  systems,  one  might  expect  a  large  frac¬ 
tion  of  PN  central  stars  to  have  gravitationally  bound  companions.  Ad¬ 
ditionally,  late-type  stars  are  more  numerous,  and  due  to  their  low  lu¬ 
minosities  would  be  difficult  to  detect  by  visual  observational  methods 
at  distances  where  PN  are  typically  found  (21  kpc).  (Only  5  known  PN 
are  thought  to  be  nearer  than  0.5  kpc).  It  is  thus  possible,  and  in  our 
view  probable,  that  a  significant  number  of  PN  nuclei  possess  cool  com¬ 
panions,  hitherto  undetected. 

Since  K  and  M  stars  emit  a  considerable  fraction  of  their  energy  in 
the  infrared,  we  calculated  the  feasibility  of  detecting  them  with  the 
Wyoming  InSb  photometer,  and  found  that  such  detections  might  be  possi¬ 
ble  out  to  distances  of  =2  kpc.  To  date  13  objects  have  been  observed, 
including  a  sample  of  planetaries  with  known  binary  nuclei  to  test  our 
hypothesis.  The  results  are  listed  below. 


Object 

IR  Excess  Remark 

Object 

IR  Excess  Remark 

NGC 

7293 

no 

NGC 

6905 

no 

NGC 

6853 

no  detection 

NGC 

6543 

yes 

known  binary 

NGC 

246 

yes 

known  binary 

NGC 

6572 

yes 

known  binary 

NGC 

7008 

? 

NGC 

6790 

yes 

NGC 

40 

7 

NGC 

6210 

yes 

new  binary? 

NGC 

6826 

yes 

new  binary  ? 

A 

63 

yes 

known  binary 

NGC 

7009 

yes 

new  binary  ? 

The  observations  were  made  using  standard  infrared  techniques  in 
the  J,  H,  K,  and  L  bands.  We  used  a  5  arc-sec  aperture  centered  on  the 
hot  star.  Any  companion  located  within  2  arc-sec  of  the  nucleus  should 
have  been  within  our  beam.  Thus,  at  a  distance  of  1  kpc,  any  binary 
system  with  a  separation  of  £  4000  A.U.  should  have  been  observed. 

We  caution  that  the  results  presented  here  are  preliminary,  and 
still  subject  to  analysis.  Specific  models  to  fit  the  data  are  present¬ 
ly  being  constructed.  We  think  that  we  may  have  detected  3  new  binary 
nuclei,  NGC  7009,  NGC  6210,  and  NGC  6826.  The  interpretation  for  NGC 
6790  is  not  clear.  This  work  was  supported  by  the  U.S.  Air  Force  Office 
of  Scientific  Research,  the  NSF,  and  University  of  Wyoming. 
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ABSTRACT 

It  is  proposed  that  cosmic  rays  may  be  the  primary  source  of 

the  incremental  energy  required  to  produce  the  bipolar  or  elongated 

structure  observed  in  many  planetary  nebulae.  If  their  energy  flux 

density  is  of  the  same  order  as  that  in  the  solar  neighborhood,  then 
4 

in  a  few  times  10  years  the  total  energy  carried  by  cosmic  rays 
through  a  region  in  space  occupied  by  an  expanding  planetary  nebula 
is  comparable  to  the  kinetic  energy  of  expansion  of  the  nebula. 
Specific  mechanisms  for  the  transfer  of  cosmic  ray  energy  to 
planetary  nebulae  are  described. 
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I.  Introduction 

Morphological  studies  of  planetary  nebulae  (PN)  indicate  that 
a  significant  fraction  of  extended  PN  exhibit  noticeable  elongation, 
or  bipolar  character.  The  long  axes  of  these  PN  have  been  shown  to 
lie  preferentially  parallel  to  the  galactic  equator  (Melnick  and 
Harwit,  1975),  and  to  the  interstellar  magnetic  field  as  determined 
from  polarization  measurements  on  stars  in  their  vicinity  (Grinin  and 
Zverva,  1968).  It  is  difficult  to  say  why  similar  studies  reported 
by  Gurzadyan  (1970)  and  Khromov  and  Kohoutek  (1968)  failed  to  show 
any  correlation  between  the  orientation  of  PN  axes  and  the  galactic 
equator,  since  few  details  of  the  investigations  were  presented  in 
either  case. 

Kirkpatrick  (1976)  has  proposed  that  the  symmetries  observed  in 
PN  are  derived  from  associated  symmetries  present  in  the  precursor 
stars,  possibly  related  to  rotation.  Livio  (1980)  thinks  that  at 
least  some  observed  symmetries  can  be  explained  by  PN  formation  around 
close  binary  central  stars.  These  mechanisms  may  very  well  operate 
in  many  instances,  but  if  the  alignment  of  PN  axes  of  symmetry  with 
the  galactic  magnetic  field  is  true,  then  at  least  one  mechanism  in 
operation  to  produce  these  symmetries  must  be  associated  with  the 
galactic  magnetic  field  in  some  way,  since  there  is  no  evidence  that 
the  axes  of  stellar  rotation  nor  of  binary  orbits  are  not  randomly 
oriented  with  respect  to  the  galaxy  (Slettebak  1949,  van  de  Kamp 
1958,  Huang  and  Wade  1966,  Huang  1967). 

Cosmic  rays  (c.r.)  are  energetic  charged  particles  (mostly 
protons)  which  naturally  follow  the  interstellar  magnetic  field 
in  their  travels,  as  they  gyrate  in  helical  paths  along  magnetic 
field  lines.  It  has  been  found  that,  through  collective  effects. 
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cosmic  rays  can  stimulate  the  formation  of  Alfven  waves  in  the  ISM 
through  which  they  travel,  and  via  this  intermediary  can  deposit 
significant  amounts  of  energy  inthe  interstellar  gas  (Parker  1968, 
1969;  Meyer  1969;  Wentzel  1972,  1974,  1977;  Ginzburg  et  al.  1973). 

In  section  II  we  show  that  if  less  than  1%  of  energy  carried 
by  cosmic  rays  is  deposited  in  a  PN,  it  is  sufficient  over  a  PN 
lifetime  to  account  for  the  observed  elongation  of  many  PN. 
Calculations  are  presented  in  section  III  for  two  possible  specific 
energy  transfer  mechanisms,  followed  by  a  brief  discussion  of  boundary 
effects  at  the  PN/ISM  interface.  We  conclude  in  section  IV  that  c.r. 
can  ultimately  be  responsible  for  the  observed  axial  symmetries  of 
at  least  some  PN. 
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II.  Energy  Considerations 


The  total  kinetic  energy  of  expansion  of  a  "typical"  PN  is 
easily  calculated.  If  we  assume  that  a  typical  nebula  contains  a  mass 
m  =  0.2  (Osterbrock  1974),  and  is  expanding  uniformly  with  a  radial 
velocity  vgx  =  20  km  sec"^  (Osterbrock  1974),  then  the  kinetic 
energy  of  expansion  is 


_  1  2  0  , a44 

T  =  o  mv  -  8  x  10  erg 
ex  2  ex  3 

Suppose  that  we  are  able  to  couple  the  c.r.  energy 

to  a  PN  with  some  efficiency  e.  Then  the  energy  transferred  from 

the  c.r.  to  the  nebula  over  some  peribd  of  time  -rcan  be  expressed  as 


/T 

E  (t)  =  e  Fc  o  (t)dt 
o 

where  Fc  =  cosmic  ray  flux 

c(t)  =  effective  nebular  cross-section  for  c.r. 
interaction 

For  a  circular  geometric  cross-section 

o(t)  =  ur2(t)  =  7r(vext  )2 

The  energy  which  may  be  transferred  during  a  period  of  time  t  is 
then 

E(t)  =  e  Fc  it  vj  t2dt 
o 


If  we  assume  that  the  coupling  efficiency  c  is  independent  of  time, 
we  obtain 


E(x)  =  !e  Fcvex  T' 


-3  -2  -1 

Using  the  c.r.  flux,  Fc  =  7  x  10  erg  cm  sec  ,  given  by  Allen  (1973), 
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4 

over  a  period  of  3  x  10  years  the  energy  transferred  to  a  PN  is 

E(x  =  3  x  104yr)  =  2.5  x  104^  e  erg 

which  is  significantly  greater  than  the  total  kinetic  energy  of 
expansion  of  the  nebula. 

Starting  with  uniform  spherically  symmetric  expansion  as  assumed 
above,  how  much  additional  kinetic  energy  is  required  to  produce  a 
PN  with  typical  elongation  (Khromov  and  Kohoutek  1968)  of  a/b  =  1.5? 
If  the  average  expansion  velocity  along  the  semi-minor  axis  b  is 
20  km  sec  \  then  the  average  expansion  velocity  along  the  semi¬ 
major  axis  a  must  be  30  km  sec  \  We.  can  create  a  crude  model  for 
the  differential  expansion  by  assuming  that  90%  of  the  nebular  mass 
continues  to  expand  at  20  km  sec'1,  and  10%  expands  at  30  km  sec'1; 
the  additional  kinetic  energy  of  expansion  is 

ATex  =  1  (0.'lm)[(vex  +  Av)2  -  ve>(2]  =  1  x  1044  erg 

Thus  the  energy  available  from  cosmic  rays  appears  to  be  sufficient 
to  account  for  the  observed  elongation  of  PN.  The  question  remains 
as  to  whether  such  coupling  can  be  shown  to  exist.  We  address  this 
problem  in  the  next  section. 


,-J.  X  . 
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III.  Coupling  Between  Cosmic  Rays  and  Planetary  Nebulae 

One  of  the  modern  discoveries  made  concerning  cosmic  rays  was 
that  of  collective  effects  (see  e.g.  the  review  by  Wentzel,  1974). 

As  c.r.  stream  through  the  galaxy  along  magnetic  field  lines,  they 
are  able  collectively  to  generate  hydromagnetic  waves  in  the  ISM 
which  propagate  through  the  gas  at  Alven  velocity 

vA  =  B  (4tt  p.)'1/2  ;  pi  =  Mini 

=  ion  mass  density 

The  Alfven  waves  are  low  frequency  transverse  electromagnetic  waves 
which  travel  along  B  parallel  to  the  c.r.  streaming  direction.  These 
waves,  in  turn,  are  capable  of  scattering  the  c.r.  which  produce  them. 

In  the  process  of  being  scattered,  the  c.r.  transfer  energy  and 
momentum  to  the  gas.  This  provides  the  basis  by  which  c.r.  energy 
can  be  transferred  to  PN.  In  the  following  we  describe  two  mechanisms 
for  c.r.  energy  transfer  to  PN,  and  comment  upon  a  third. 

A.  Direct  Interaction 

Following  Wentzel  (1974),  consider  cosmic  rays  streaming  down 
a  c.r.  density  gradient  through  an  ionized  region  (PN).  We  have 

4 

seen  in  section  II  that  the  energy  available  from  c.r.  over  t  =  3  x  10  yr 

is  2.5  x  lO^eerg,  and  the  energy  required  for  elongation  of  PN  is 
44 

10  erg.  The  efficiency  of  the  energy  transfer  mechanism  must 
then  be  e  =  .004.  The  c.r.  energy  loss  over  a  distance  of  one  c.r. 
density  scale  height  (L)  is  given  as 


where  <v>  is  the  c.r.  streaming  velocity  in  the  medium.  For  intercloud 
space  Wentzel  (1974)  gives 
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<v>  =  55  +  nO(-n—  )1,5  km  sec"1 


so  that  for  c.r.  momentum  p-  Mc  (~1  GeV  c.r.),  <v>  =165  km  sec'1. 
In  a  highly  ionized  region  (e.g.  PN),  <v>  is  given  by  Wentzel  as 


1  ,,  ,100  pc  ,'72  ,  P  ,0'75 

<v>  =  3  Y  Vfl  +  56  (— [K-)  (-^) 


km  sec 


Using  typical  values  suggested  by  Wentzel  for  the  c.r.  spectrum  index 
Y  =  4.5,  the  c.r.  density  scale  height  L  =  100  pc,  and  p=  M  ,  we 
find  <v>  =1.5v^  +  56  km  sec'1,  placing  a  lower  limit  on  the  streaming 
velocity  of  56  km  sec'1  in  a  PN,  with  the  exact  value  depending  on  the 
magnetic  field  strength  B  in  the  ionized  region,  through  the  Alfven 

O  _  O 

velocity  v^.  The  Alfven  velocity  inside  a  typical  PN  (ni  =10  cm  ) 
will  be  vft  =0.2  km  sec'1  if  the  magnetic  field  strength  is  the  same 
as  the  interstellar  value  (B  =3  x  10  ^  G),  or  vft  =  20  km  sec'1  if 
B  should  happen  to  be  2  orders  of  magnitude  larger. 

The  time  required  for  c.r.  to  stream  down  one  c.r.  density 
scale  height  L  is  t  =L/<v>  *  10^  yr,  within  a  factor  of  2,  depending 
on  <v>.  During  the  time  t  =10^  yr  c.r.  will  transfer  AE  =  (vA/<v>)E  erg 
to  the  medium  through  which  they  pass,  and  therefore,  the  energy 

4 

transferred  during  a  PN  lifetime  (x  =  3  x  10  yr)  is 


3  x  10^ 


E  erg 


=  °-03  TV>  E  erg 


In  order  for  this  mechanism  to  be  viable  we  require 


e  «  .004  a  .03 
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Inside  the  ionized  region  we  showed  that  <v>  ~1.5vA  +  56  km  sec"1. 

Combining  this  with  the  above,  the  Alfven  velocity  in  the  PN  must 

be  i 

vA  -  9  km  sec 

The  corresponding  magnetic  field  which  must  exist  in  the  PN  is  then 
B  *  10"4  Gauss 

in  order  for  substantial  c.r.  energy  to  be  deposited  in  the  nebula. 

We  have  thus  far  implicitly  assumed  that  the  energy  of  cosmic  ray 

heating  is  coverted  to  kinetic  energy  of  expansion.  Such  is  not 

necessarily  the  case,  however,  as  we  now  show.  For  a  typical  electron 

kinetic  temperature  of  T  s  1  x  104  K‘,  the  speed  of  sound  in  a  PN  will 

be  cQ  =  (2NgkTe/p)^  -  13  km  sec’1.  The  theory  of  expansion  of  a  gas 

cloud  into  vacuum  as  a  result  of  gas  pressure  is  discussed  by  Osterbrock 

(1974),  the  main  result  applicable  here  being  that  the  bulk  of  the  gas 

expands  outward  at  a  velocity  slightly  greater  than  the  sound  velocity, 

cQ.  Thus,  in  order  for  expansion  velocities  on  the  order  of  30  km 

sec’1  to  be  attained,  the  cosmic  rays  must  heat  the  nebular  gas  to 
2  4 

Tg  =  pcq  /2Ngk  -  5  x  10  K.  The  radiative  cooling  rate  in  a  typical  PN 

(T  ~  104  K)  is~3x  10’24  Ng2  erg/cm^  -  sec  (Osterbrock,  1974).  If  ~10& 

of  the  nebular  gas  is  heated  by  c.r.  as  discussed  above,  then  the  c.r. 

-24  2  3 

heating  rate  will  be  ~1  x  10  N  e  erg/cm  -sec.  The  additional 
energy  imput  by  c.r.  will  result  in  increased  heating  and  cooling  rates, 
with  a  maximum  increase  in  the  equilibrium  kinetic  temperature  of  ~900  K, 
for  the  case  ofe=  1.  That  is  clearly  insufficient  to  produce  the 
required  expansion  velocities  of  '-30  km  sec"1.  It  therefore  appears  that 
most  of  the  c.r.  energy  deposited  in  a  PN  by  direct  interaction  will  be 
converted  to  photons  and  radiated  away  rather  than  going  into  expansion. 
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B.  Interstellar  Medium  as  an  Intermediary 

We  now  consider  a  process  whereby  the  ISM  is  heated  by  c.r.  as 

they  stream  through  it.  The  hot  ionized  portion  of  the  ISM  is  then 

constrained  to  move  along  galactic  magnetic  field  lines  and  upon 

encounter  with  a  PN  will  preferentially  heat  those  regions  of 

the  PN  into  which  the  galactic  magnetic  field  lines  lead,  thus 

depositing  the  requisite  kinetic  energy  in  the  PN  at  the  right  places 

to  explain  the  observed  elongations. 

Savage  and  de  Boer  (1979)  have  reported  observational  evidence 

for  the  existence  of  a  hot  gaseous  galactic  corona  which  they  con- 

.  5 

elude  most  likely  has  a  kinetic  temperature  of  ~10  K.  We  will 
now  show  that  this  can  be  explained  in  terms  of  c.r.  heating. 

As  c.r.  stream  through  the  ISM  some  scattering  occurs,  and 
energy  and  momentum  are  deposited  in  the  ISM.  The  scattering 
efficiency  will  depend  on  the  c.r.  streaming  velocity  <v>  and  the 
Alfven  velocity  v^  in  the  medium.  For  intercloud  space  <v>  =165  km  sec 
(sec.  II  A.).  The  Alfve'n  velocity  in  the  ISM  is  v^  =38  km  sec 
for  B  =  3  x  10  and  n..  =3  x  10  ^  cm  ^ (Gomez-Gonza'lez  and  Guelin 
1974).  The  enerq.v  transferred  to  the  ISM  during  the  time  to  traverse 
one  scale  height  (t=  L/<v>=6  x  10^  yr)  is  then 

VA 

AE  =  E=  0.23E 
<  v  > 


Suppose  the  c.r.  flux  is  more  or  less  constant  over  the  time  t  in  a 


given  region  of  space;  then  the  energy  per  unit  volume  deposited 

5 

by  c.r.  in  the  ISM  in  6  x  10  yr  will  be 

0.23F 


AE 

V 


<  v  > 


c  1n-10  -3 

=10  erg  cm 


or,  assuming  equipartition  of  energy  between  electrons  and  ions, 
the  energy  transferred  to  each  particle  of  the  ISM  is 
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Now  suppose  that  at  some  point  in  time  a  PN  starts  to  form 

5 

in  a  region  of  space  which  has  been  heated  to  '-10  K  by  cosmic 
rays.  We  note  that  the  hot  ISM  ions  will  be  constrained  to  move  along 

5 

interstellar  magnetic  field  lines,  since  at  T  =10  K  the  gas  pressure 

? 

nkT  is  expected  to  be  comparable  to  the  magnetic  pressure  B  /8ir, 
rather  than  overwhelming  it.  The  hot  ISM  ions  will  therefore  approach 
a  PN  primarily  along  magnetic  field  lines. 

Some  net  momentum  associated  with  the  c.r.  streaming  will  have 
been  transferred  to  the  ISM,  and  subsequently  to  any  PN  encountered 
by  the  ISM  ions.  However,  we  must  realize  that  the  net  momentum  is 
the  vector  sum  of  momenta  of  nearly  equal  numbers  of  ions  travelling 
in  opposite  directions  along  the  magnetic  field.  The  ISM  can  there¬ 
fore  interact  with  and  deposit  energy  simultaneously  in  opposite 
sides  of  a  PN.  The  deposited  energy  results  in  local  heating  and 
expansion  of  the  nebular  gas,  producing  the  characteristic  axial 
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symmetries  observed  in  many  PN. 

Now  consider  the  interaction  of  the  ions  of  the  ISM  with 
those  of  a  much  more  dense  PN  which  they  encounter.  This  will 
be  characterized  by  Coulomb  collisions,  and  the  ion-ion  interactions 
can  be  described  in  terms  of  an  effective  momentum  transfer  cross- 
section  (see  Appendix  A) 


2  An  m.v.‘ 

0 .  - «" J c-Vt2- 

C  7  P 

m.v.  e 

i  i 


where  Z  =  charge  no.  =  1  for  hydrogen  ions 
e  =  4.8  x  10  ^  esu 
m.  =  1.67  x  10'24  g 
v.  =  50  x  10^  cm  sec’^ 

Aq  =  22  cm  =  Debye  shielding  length  in  PN  (Appendix  A) 

-15  2 

We  find  am  =8.5  x  10  cm  .  The  corresponding  mean  free  path 
(mfp)  for  the  ISM  ions  in  the  PN  is 

mfp  = — —  =10^  cm  =4  x  10  ®  pc 

n.o 
l  m 

where  in  this  case  n..  =  ion  density  in  PN. 

Another  approach  to  exmaining  this  interaction  is  to  calculate 
the  "slowing  down"  time,  t$,  given  by  Spitzer  (1956). 

4 

We  find  t  -1.3  x  10  sec,  with  the  corresponding  mfp  for  50  km 
sec’l  ions 

mfp  =  v^.t^  =  2  x  10'^  pc 

in  good  agreement  with  the  momentum  cross-section  calculation. 
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We  see  that  the  mfp  for  incident  ions  of  the  ISM  is  small  with 

respect  to  the  dimensions  of  the  nebula  at  nearly  every  stage  of 

expansion.  Therefore  the  energy  carried  by  the  ISM  must  be  completely 

r  transferred  to  the  ions  of  the  PN.  It  is  important  to  realize  that 

the  c.r.  energy  thus  depositied  will  result  in  purely  local  heating. 

-8  “6 

For  a  mfp  of  ~3  x  10  pc  calculated  above,  only  ~10  of  the  nebular 

mass  will  receive  c.r.  energy.  The  c.r.  heating  rate  in  the  local 

-19  2  3 

regions  of  interaction  will  then  be  ~1  x  10  Ng  e  erg/cm  -sec.  The 

kinetic  energy  of  a  50  km  sec”*  hydrogen  ion  is  2  x  10"11  erg,  which 

gives  an  efficiency  factor  t  =  (0.23)12  x  10”**/1.6  x  10”^)  =  .003  for 

transfer  of  c.r.  energy  to  the  ISM  and  subsequently  to  PN.  The  c.r. 

-22  2  3 

local  heating  rate  in  PN  is  then  ~3  x  10  Ng  erg/cm  -sec,  a  factor 

4 

of  -100  greater  than  the  radiative  cooling  rate  for  a  10  K  nebula. 

Most  of  the  c.r.  energy  will  therefore  be  available  to  be  converted  to 
kinetic  energy  of  expansion,  with  very  little  of  it  lost  to  radiative 
cooling.  The  c.r.  energy  thus  transferred  during  a  typical  PN  lifetime 
(3  x  104  yr)  will  then  be 

E(t  =  3  x  104  yr)  =  2.5  x  1046e  erg  (sec.  II) 

=  7.5  x  1043  erg 

Even  though  the  calculations  have  been  quite  crude,  we  see  that  this 

44 

result  indicates  reasonable  agreement  with  the  10  erg  of  additional 
energy  required  to  explain  PN  elongation  (sec.  II). 

The  detailed  scenario  of  PN  interaction  with  the  ISM  over  a  PN 
lifetime  is  most  likely  quite  complicated,  and  is  not  intended  to  be 
the  subject  of  this  paper.  We  note  that  as  the  outer  layers  of  the 
PN  expand  as  a  result  of  the  c.r.  heating,  the  electron  density  Ng 
decreases,  and  the  mfp  therefore  increases.  In  this  way  deeper  layers 
of  the  PN  are  able  to  interact  with  the  ISM  ions,  and  it  is  reasonable 
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to  expect  that  eventually  a  significant  fraction  of  the  PN  mass  will 
be  subject  to  c.r.  heating,  the  energy  of  which  will  be  available  for 
subsequent  expansion. 

C.  Boundary  Effects 

It  may  be  that  the  sharp  boundaries  which  exist  between  the 
ISM  and  PN  can  contribute  significantly  to  the  energy  transfer. 

The  configuration  of  the  interstellar  magnetic  field  in  the 
vicinity  of  PN  is  not  known.  However,  it  is  well-known  that  the 
diffusion  of  a  magnetic  field  through  astrophysical  plasmas,  such 
as  those  we  are  considering  here,  is  extremely  slow,  so  much  so 
as  to  effectively  "freeze"  the  field  into  the  plasma.  The  expanding 
motion  of  the  PN  plasma  would  presumably  then  carry  any  embedded 
magnetic  field  with  it,  since  the  gas  pressure  is  several  orders 
of  magnitude  larger  than  the  magnetic  pressure.  This  would  produce 
some  sort  of  kinking  or  bending  of  the  magnetic  field  at  the  PN/ISM 
boundary.  Theoretically,  fluctuations  in  the  galactic  magnetic 
field  with  dimensions  on  the  order  of  10  pc  will  scatter  cosmic 
rays  (Wentzel  1974),  thereby  transferring  energy  from  the  c.r.  to 
the  surrounding  medium.  Whether  variations  in  the  magnetic  field  of 
this  magnitude  are  present  near  PN/ISM  boundaries  is  not  known,  but 
the  possibility  exists  and  could  be  important  for  transferring  c.r. 
energy  to  PN. 
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IV.  Summary 

Attempts  to  invoke  magnetic  fields  (Gurzadyan  1962,  1970; 

Woyk  1968)  in  direct  interaction  with  planetary  nebulae  have  not 
been  fruitful  (Kahn  1968;  Wolt.ier  1968;  Menzel  1968).  Castinq 
the  galactic  magnetic  field  in  an  accessory  role,  as  proposed  here, 
seems  to  be  more  promising.  Cosmic  rays  of  1  GeV  or  more  have  been 
shown  to  carry  sufficient  energy  to  supply  an  expanding  PN  with 
the  additional  energy  required  to  produce  the  observed  elongat’  ,ns. 
Collective  effects  from  cosmic  ray  theory  have  been  used  to  explain 
how  energy  may  be  transferred  from  the  cosmic  rays  to  the  PN.  It 
seems  that  sufficient  energy  to  explain  the  observations  may  be  gained 
by  the  PN  from  interaction  with  the  interstellar  medium  which  has  been 
heated  by  c.r.,  although  other  possibilities  are  not  ruled  out. 

The  role  which  boundary  effects  might  play  remains  to  be  considered 
in  more  detail.  The  existence  of  magnetic  fields  within  PN  of 

sufficient  strength  to  make  direct  interaction  of  c.r.  with  nebular 
material  a  viable  energy  transfer  mechanism  remains  yet  to  be 

disproved.  Observations  of  Zeeman  broadening  in  high-n  hydrogen 
(radio  re combination)  lines  might  be  used  (Greve  1975,  Troland 
and  Heiles  1977,  Greve  and  Pauls  1980)  to  detect  or  set  an  upper 
limit  to  magnetic  field  strength  in  PN. 

Cosmic  rays  streaming  along  magnetic  field  lines  may  be  the 
primary  source  of  incremental  energy  for  elongation  of  expanding 
planetary  nebulae.  This  is  consistent  with  our  present  understanding 
of  cosmic  ray  physics,  observations  of  the  hot  gaseous  galactic 
corona,  and  alignment  of  planetary  nebula  axes  along  the  galactic 
magnetic  field. 
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The  fractional  change  in  momentum  along  the  forward  direction  will 
be  proportional  to  (1-cosX).  This  leads  us  to  define  a  momentum 

transfer  cross-section,  om,  as  follows 

'  v 

am  =  (1  -  cos  x)  2tt  sin  xdx 

•  x  *  *  min 

Using  the  Rutherford  formula,  this  becomes 


By  trigonometric  manipulation  this  can  be  put  in  the  form 
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,  r  -|2 

=  4tt  [ - JT’J 


Vi 


— ■  -X/-—  d  (x/2) 


sin  x/2 
min 


The  solution  is  obtained  by  letting  y  -  sin  x/2,  and  the  result  is 

2 

-]  In  A 


qiqs  -2 

a  =  4v  [ - 

m 


Vi 


where  A  = 


1 


/  X  minx 
sin  (  o  / 


We  note  that 


1  im  A  = 
X  -*■  o 


X  min 
2 


b  2 

max  m. v . 

i  l 


qiqs 


In  a  plasma  it  seems  reasonable  to  choose  the  maximum  impact 
parameter,  b  .  ,  equal  to  the  Debye  shielding  length,  An 

fTlaX  U 


max 


4Tinee 


In  a  PN  for  which  Te  =10^K  and  ng  =-703  cm"3,  we  find  A  ^  =  22cm. 
Letting  q.  =  qg  =  Ze,  the  momentum  transfer  cross-section  takes 
the  final  form 


2  2 

a  =  47T  [-!§ - ] 

m  L  o 

mi  vi 


A0  m.v. 


'Ze 


2 
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